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Abstract. 

The current status of tests of the theory of strong interactions, Quantum Chromo 
Dynamics (QCD), with data from hadron production in e"'"e~ annihilation experiments 
is reviewed. The LEP experiments ALEPH, DELPHI, L3 and OPAL have pubhshed 
many analyses with data recorded on the T? resonance at ^fs = 91.2 GeV and above 
up to y/s > 200 GeV. There are also results from SLD at \/s = 91.2 GeV and from 
reanalysis of data recorded by the JADE experiment at 14 < < 44 GeV. The 
results of studies of jet and event shape observables, of particle production and of quark 
gluon jet differences are compared with predictions by perturbative QCD calculations. 
Determinations of the strong coupling constant as('7iz°) from jet and event shape 
observables, scaling violation and fragmentation functions, inclusive observables from 
TP decays, hadronic r decays and hadron production in low energy e"'"e~ annihilation 
are discussed. Updates of the measurements are performed where new data or improved 
calculations have become available. The best value of 03(^120) is obtained from an 
average of measurements using inclusive observables calculated in NNLO QCD: 

ag(mzo) = 0.1211 ±0.0021 
where the error is dominated by theoretical systematic uncertainties. The other 
measurements of ag(m2o) are in good agreement with this value. Finally, investigations 
of the gauge structure of QCD are summarised and the best values for the colour factors 
are determined: 

Ca ^ 2.89 ±0.21 
Cf = 1.30 ±0.09 

with errors dominated by systematic uncertainties and in good agreement with the 
expectation from QCD with the SU(3) gauge symmetry. 
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1. Introduction 

The interactions between the constituents of matter are successfully described by the 
four forces: the weak, electromagnetic and strong forces and the gravitational force. 
The weak and the strong interactions occur at small atomic to subatomic distances, 
the electromagnetic interaction is observed at subatomic to macroscopic distances while 
effects of Gravitation only play a role at macroscopic distances. 

The strong interaction, the main focus of this review, is responsible for the 
existence of all composite elementary particles (hadrons) by providing the binding 
force between the constituents and also for most of the short lived hadron decays. 
Furthermore, the binding of protons and neutrons in nuclei may be explained in analogy 
to chemical binding of molecules based on the strong interaction of the proton and 
neutron constituents. 

The constituents of hadrons are known as partons or quarks. The known spectrum 
of hadrons can be explained as composites of five (out of a total of six) quark flavours, 
where each flavour occurs in three variants distinguished by the so-called colour. In 
tabled the basic properties of the six known quarks are shown. 

Table 1. Basic properties of the six quarks in the standard model 



Generation 1st 


2nd 


3rd 


Charge 


up-type u (up) 


c (charm) 


t (top) 


-h2/3 


down-type d (down) 


s (strange) 


b (bottom) 


-1/3 



A dynamic theory of strong interactions at the constituent level. Quantum Chromo 
Dynamics (QCD), is constructed as a renormalised field theory in close analogy to 
Quantum Electro Dynamics (QED), the quantum field theory of the electromagnetic 
interaction, see e.g. [1-6]. The quarks are viewed as carriers of three different strong 
charges, referred to as colours. Interactions between colour charged quarks are mediated 
by gluons, in analogy to the exchange of photons in QED. The theory is constructed 
to remain invariant under exchange of the three colour charges, i.e. local gauge 
transformations in the colour space with the SU(3) symmetry. From the requirement 
of local gauge invariance under SU(3) the special properties of the gluons are derived: 
there are eight gluons each carrying a colour charge and anti-charge and thus the gluons 
can directly interact with each other. The gluon-gluon interactions of QCD are not 
present in QED; QCD is referred to as a non-abelian gauge theory while QED is an 
example of an abelian gauge theory. For massless quarks the strong coupling constant 
ag is the only free parameter of the theory. 

Qualitatively, some important properties of QCD already follow from the possibility 
of direct gluon-gluon interactions. In QED, the scale dependence of the coupling 
constant a, the so-called running, may be seen as a consequence of screening of the bare 
electric charge due to vacuum polarisation. At small momentum scales, the exchanged 
photon resolves only a large volume around the bare charge and thus is exposed to 
stronger charge screening. At large momentum scales a smaller volume is resolved 
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and less charge screening occurs leading to a rising value of the coupling constant. In 
QCD the effect has the opposite direction, because a cloud of virtual colour charged 
gluons around a bare colour charge effectively spreads the colour charge over a volume 
around the quark carrying the bare colour charge. Thus at small momentum scales 
exchanged gluons interact with an apparently stronger colour charge while at large 
momentum scales a correspondingly smaller fraction of the colour charge is resolved. 
In the limit of infinite momentum scales the resolved charge would vanish, leading to 
the prediction of asymptotic freedom for quarks. At very small momentum scales the 
resolved colour charge becomes so large that quarks are tightly bound into hadrons; this 
explains the observation that quarks are confined into hadrons and cannot be observed 
as free particles in an experiment. 

The process of electron positron annihilation into hadrons, e"'"e~ — >■ hadrons, is an 
ideal laboratory for studies of strong interaction phenomena at the parton, i.e. quark 
and gluon level, for several reasons: 

• There is no interference between the initial and the final state with the consequence 
that interpretation of the hadronic final state in terms of QCD processes is 
simplified. 

• The four-momentum of the initial state is, in the absence of significant QED initial 
state radiation (ISR), fully transferred to the final state. In most experiments the 
electron and positron beams have equal energies such that the centre of mass system 
(cms) of the final state coincides with the laboratory frame. 

• The experimental conditions in e"'"e~ annihilation are usually rather clean, because 
in past and existing colliders there were no multiple interactions in a given bunch 
crossing and backgrounds coming directly from the accelerator like lost beam 
particles or synchrotron radiation were kept at a low level in the experiments. 

• The presently available e"'"e~ annihilation data allow detailed studies of QCD in a 
large range of cms energies covering more than an order of magnitude from 12 GeV 
to 209 GeV. This unique collection of data forms in particular the basis of tests of 
the energy scale dependence of QCD predictions. 

Figure H presents an overview of our current understanding of the process e"'"e~ — *• 
hadrons. The incoming electron and electron annihilate into an intermediate vector 
boson, shown as (Z°/7)*, possibly after initial state radiation. The quark pair produced 
in the vector boson decay will start to radiate gluons which in turn radiate gluons or 
turn into quark pairs themselves. In this way a so-called parton shower develops until 
all parton interactions happen at low energy scales of about 1 GeV where confinement 
and thus hadron formation sets in. The parton shower occupies a volume with radius 
~ 10"^ fm, according to the uncertainty principle. 

This review will begin in section |2l with a brief introduction into fundamentals of 
QCD predictions and sect ion IHl gives a short overview of milestones in the history of QCD 
tests in e"'"e~ annihilation together with a summary of the e"'"e~ colliders and experiments 
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Electro-weak Production Parton Shower Hadronisation 




Vs' ~ 1 GeV 



Parton Level Hadron Level 

Figure 1. Schematic view of a e+e" liadrons event; \/s' denotes the invariant mass 
of the hadronic system, usually close the nominal cms energy ^/s of the experiment. 

which produced the currently available data. The following section HI discusses studies 
based on jet reconstruction algorithms and event shape observables. The subject of 
section El is inclusive particle production and section IHl presents results from inclusive 
observables. A summary of determinations of the strong coupling as is given in sectional 
while the properties of jets originating from quarks or gluons are discussed in section |H1 
Studies of the gauge structure of QCD are reviewed in section IHl followed by conclusions 
and outlook in section ITUl 

2. Basics of QCD and hadron production 

We collect here the important predictions of QCD which can be studied with data from 
hadron production in e"'"e~ annihilation. In QCD with the SU(3) gauge symmetry the 
following fundamental processes are possible [5,6]: 

Gluon radiation from quarks This process is the analog of photon bremsstrahlung 
and occurs with relative strength Cp = 4/3. 

Quark pair production This process is the analog of e.g. e"'"e~ pair production from a 
photon and has a relative strength Tp = 1/2. However, it is possible to create quark 
pairs of all kinematically accessible quark flavours Uf and thus the contribution of 
this process is effectively TpUf. 

Triple gluon vertex This process is unique to QCD because the gluons are themselves 
colour charged, the relative strength is Ca = 3. 
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There is in addition the quartic gluon vertex which is at least 0{a'^ in e+e~ annihilation 
and known to be a negligible contribution to radiative corrections to lower order 
processes [7]. Figure El shows the Feynman diagrams of the three fundamental processes. 




Figure 2. Feynman diagrams of the tliree fundamental QCD processes. 

The overall strength of the strong coupling is given by a constant of nature, the 
strong coupling constant as. This constant is the only free parameter of the theory when 
quarks are treated as massless. Therefore many studies of QCD can be formulated as 
measurements of the strong coupling constant and comparisons of the results for as of 
different analyses give insight into the consistency of the theory. 



2.1. Everything runs 

2.1.1. Renormalisation scale dependence The predictions of QCD are finite, because 
the theory has been renormalised in order to remove divergent terms. The 
renormalisation introduces a dependence of QCD predictions on the energy scale /z, 
at which the renormalisation is performed. The renormalisation scale /i is arbitrary and 
physical results should be independent of /i. Several prescriptions for renormalisation 
exist, the so-called renormalisation schemes, and in this report we use the widely adopted 
MS renormalisation scheme [5,6] throughout. 

A QCD prediction of an observable quantity R measured in e"'"e~ annihilation, e.g. 
a differential cross section, can be written as: 

oo 

R=Y.Rnal , (1) 

ra=0 

where the Rn are the coefficients in nth order of the perturbation series. The quantity 
R depends on as and the ratio x^ = Q/ oi the physical scale Q of the process and the 
renormalisation scale /i, for reasons of consistency of the dimension of R. In the limit 
n — i> oo this dependence on the renormalisation scale is expected to vanish; however, for 
truncated perturbative calculations a residual dependence will remain. The requirement 
of renormalisation scale independence is formally expressed as the renormalisation group 
equation (RGE), also known as Callan-Symanzik equation [8,9], see e.g. [5]: 
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The renormalisation scale dependence of the coefficients i?„ will be compensated by a 
renormalisation scale dependence of the coupling constant as(/i): 

Equation Q introduces the /5-function of QCD, with coefficients [5] 

^" 127r ' ^ ^ 

_ 17CI - 5C^n^ - SCFUf 



P2 



2857C3 + (54C^ - Q15CfCa - 1415Cl)n/ + (66Cp + 79CaW 



f 




34567r3 

The coefficient (32 depends on the renormalisation scheme and is shown here for the 
MS scheme. The variable Uf specifies the number of active quark flavours which are 
considered in the calculation and thus from kinematics it depends on the energy scale 
of the process. At low energies below the thresholds for heavy quark production direct 
contributions as well as contributions to radiative corrections are suppressed resulting 
in Hf = 3. When with increasing energy scale a heavy quark production threshold is 
crossed nj must be increased by one unit, see section !^. 1.41 below. 
Combining equations ((H), ^ and (jH)) results in 

= -"'S + + 4if^') (-"'P - RiPo^ + (5) 

A solution for equation (jSJ can be found by demanding that the coefficients of ag vanish 
for all orders n. After integration from a specific choice of renormalisation scale to 
another scale the coefficients can be written as 

Rq = const. , (6) 

-Ri = const. , 

R2ixl) = R2{1) + PoRi Inxl and 

RsixD = Rs{l) + (2i?2(l)/3o + RiPi) InxJ + R^P^, In^ 

with = Q/fiR. The Ri{l) are the coefficients evaluated at the renormalisation scale 
/iR. We note that all coefficients Rn, n > 2 depend explicitly on the renormalisation 
scale parameter x^. However, the coefficients Rn{x^) only depend on x^ via terms 
containing Inx^. This observation has the important consequence that QCD predictions 
for physical observables at any scale given by can be computed in terms of constant 
coefficients derived using the renormalisation scale given at x^ = 1. The scale 
dependence of the prediction then enters only through the scale dependent strong 
coupling constant as(Q^). 

The choice of the renormalisation scheme (RS) is arbitrary and usually dictated by 
convenience to perform a calculation. It is also evident that a result of a QCD prediction 
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for a physical observable like a cross section should in principle not depend on the choice 
of RS. In practice with perturbative predictions truncated at e.g. 2nd or 3rd order there 
are significant dependences on the RS which may be interpreted as indications of the 
size of missing higher order terms. 

This problem and attempts to provide solutions have already been discussed more 
than 20 years ago when the first next-to-leading (NLO) QCD calculations became 
available [10-14]. Three different ways to choose an "optimal" renormalisation scheme 
or scale for sufficiently inclusive observables depending on a single energy scale were 
devised, see e.g. [15] for a review. Examples for such observables are mean values 
(or higher moments) of event shape distributions or -Rg+e-- It is important to note 
that at NLO a variation of the renormalisation scale is equivalent to a variation of 
the RS yielding the same results [16]. The three methods for choosing an optimal 
renormalisation scheme are the following: 

Principle of minimum sensitivity (PMS) The principle of minimum sensitivity is 
based on the observation that at a given order all possible RSs can be labelled by 
their renormalisation point (scale) and the coefficients of their /5-functions [12]. 
Now one may find the RS for which a QCD prediction -R(as) for a physical 
observable has minimal sensitivity w.r.t. choosing a RS. This corresponds to a 
stationary point of the function -R(as, h) with the RS labels /j. At NLO finding a 
stationary point of R{a^,x^) is sufficient, since at this order the coefficients of the 
/3-function are universal. 

Method of effective charges (ECH) Using the freedom to choose the RS the 
perturbative expansion is rearranged such that higher order terms vanish. The 
running coupling absorbs all scale dependent effects and becomes a process 
dependent so-called effective charge which satisfies a generalised /3-function [10,11]. 
This method is also referred to as "fastest apparent convergence" or FAC. At NLO 
a convenient way to find the FAC scheme is to find the value of for which the 
complete 0{a'^ term vanishes. 

Brodsky-Lepage-Mackenzie method (BLM) The BLM method proposes to ad- 
just the renormalisation point (scale) for a given RS such that the dependence 
of the NLO term on the number of quark flavours ^/ is cancelled [13, 17]. This 
prescription implies that vacuum polarisation corrections due to fermion pairs are 
absorbed by the running coupling as(/i). 

2.1.2. The running strong coupling constant The scale dependence or running of 
the strong coupling constant can be derived by integrating equation Q to obtain an 
expression relating values of as at two different scales Q and /xr. Since the choice of 
renormalisation scale is arbitrary the expressions can be used to relate values of as at 
any two scales. Solutions differ according to how many orders (often referred to as the 
number of loops) in perturbation theory have been considered in the calculation of the 
/3-function. The first term corresponds to one order (loop), the second to two orders 
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(loops), etc. The solutions valid for one and two loops are: 

«s(Q^) = , and (7) 



as(<5^) tts(/iR) /5o V"s(/^r) /?o + Aas((5^)/ 

The 2-loop solution shown in equation (jH)) may be solved numerically for aQ,{Q'^). The 
solution for a 3-loop /3-function is also found directly by integrating equation Q: 

/3olnxJ = F(as(Q))-i^(as(/^R)) • (9) 

The function -F(as) = Po I ~l/(/^o«i + Ao^s + /^2«s) '^'^s is for < 5, i.e. P2 > 0, given 
byf 

1 f3! - 2/30/^2 ^ ( (3i + 2(32as , 
(as) = arctan — ^^^^^^= (10) 

«s /5ov'-/5? + W2 V\/-/5? + W2 



A log(as) _ Pi log(/?o + Pias + /g2«s) 

As before equation Q may be solved numerically for as{Q)- For the case n/ > 6 we 
have < and the second term of the RHS of equation (fTUI) is rewritten to contain 
only real terms with z = {Pi + 2f32a^)/ \J Pi — '^PoP2'- 

Pf-2PoP2 1, 1 + z 
. = - log . (11) 

Po^/p[^^mh'^ 

The case Uf = 6 becomes relevant when an evolution of as to energy scales above the 
threshold for t quark production is performed. 

The perturbative running of the strong coupling constant breaks down at sufficiently 
small scales Q such that e.g. in LO as(/iR)/?o In = — 1. At such low scales values of the 
strong coupling are as > 1 and thus perturbative expansions in as will fail to converge. 
The value of Q = Aqcd where this breakdown of perturbative QCD occurs is known as 
the Landau pole. Since typical light hadron masses are (9(100) MeV one might expect 
that Aqcd has similar values. Using the LO equation ((Zj) this yields as(mzo) ^0.1 and 
as(l GeV^) ~ 0.2 which sets a rough lower limit for the apphcability of perturbative 
QCD. 



2.1.3. Running quark masses QCD predictions discussed so far are valid for massless 
quarks. However, it is well known that quarks have masses ranging from 0{1) MeV 
for the light u and d quarks to about 174 GeV for the t quark [18]. Mass effects are 
incorporated into the theory via mass terms in the Lagrangian which are subject to 
renormalisation like the coupling constant. The requirement of independence of QCD 

f Integration courtesy of integrals.wolfrain.com. 
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predictions from the renormalisation scale introduces in the presence of a quark mass a 
new term into the RGE equation (j2)): 

^ (s/i^ ~'~ dii^ das ~'~ dfi'^ drnj ^ ^ ^^'^^ 

with 

^^1^ = -7m(/^^)"^(/^^) (13) 

where ^mif^"^) is the so-called mass anomalous dimension. The mass anomalous 
dimension 7^ has been calculated as a power series in the strong coupling constant, 
7m(tts) = 7oQ^s +7iQ^s + ^("^s)' '^it^ coefficients 70 = 1/vr and 71 = (303 — 10nj)/(727r^) 
in the MS renormalisation scheme [19]. 

Integrating equation (fT^ and substituting as(/i|^) for /ip^ gives an equation for the 
running quark mass: 

^(Q')7m(«s). 

Solutions to equation |THl in the MS renormalisation scheme are e.g. given in [19]. 



m(g ) = m(/iR)exp / ———dag • (14) 



2.1.4. Heavy quark thresholds The evolution of the strong coupling constant, 
equation Q, and the quark masses, equation (dH), depends on the number of active 
quark flavours Uf through the coefficients of the QCD /3-function and mass anomalous 
dimension, respectively. When an evolution of e.g. as(/^^) across the excitation threshold 
/iq for production of heavy quark pairs is attempted an explicit treatment of the changing 
number of flavours, rif —>■ Uf + 1, is necessary. 

One requires that the theory at scales below the heavy quark threshold fiq with n/ 
active quark flavours is consistent with the theory at scales above /iq with Uf + 1 active 
quark flavours [20]. This results in matching conditions for the strong coupling constant 
at the heavy quark threshold; these are asi^f) = as{nf + 1) in LO (leading order) and 
NLO. At the next order the matching involves a discontinuity of the value of the strong 
coupling constant at the heavy quark threshold /iq: 

as(^/,/iq) = as(^/ + 1,/Uq) + ai(?i/ + l,/iq)c2 , (15) 

where the value of the coefficient C2 depends on the definition of the heavy quark mass. 
For the running MS mass, i.e. /iq = '"^qMs("^qMs)' '^2 = ll/(727r^) while for the 

pole mass Mq with /iq = Mq one has C2 = — 7/(247r^) [20]. 

In this report we use 3-loop running of the strong coupling constant as explained 
in section 1^.1.21 together with matching at heavy quark thresholds using equation ()15j] 
and MS masses '^c,ms(^c,ms) ~ 1-25 GeV and '^b,Ms("^f) Ms) ~ '^•^S GeV [18] when a 
value of as is evolved from an energy scale Qi to another scale Q2- 
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2.2. Perturhative QCD 

Predictions of perturbative QCD are possible for observables which fulfil the Sterman- 
Weinberg criteria of infrared and coUinear safety [21]: an observable is infrared and 
collinear safe when its value is not affected by the emission of low-momentum partons 
or by the replacement of a parton by collinear partons with the same total 4-momentum. 
For such observables predictions as power series in the strong coupling ag may be 
performed. 

The observables typically used in studies of e"^e~ annihilation to hadrons may be 
classified as follows: 

exclusive Exclusive observables derived from jet clustering algorithms or based on 
event shape definitions (see section |3] for details) classify the hadronic final states 
according to their topology. Measurements of differential cross sections for such 
observables may be compared with corresponding QCD predictions. 

semi-inclusive With semi-inclusive observables all hadronic final states are considered 
and average properties of the produced hadrons such as particle multiplicity or 
momentum spectra are studied. 

inclusive The inclusive observables are simply based on counting of hadronic final 
states. Examples are the hadronic widths of the Z° boson or the r lepton or the 
total cross section for hadron production in e"'"e~ annihilation. 

Most exclusive observables used in hadron production from e"'"e~ annihilation are 
sensitive to the non-collinear emission of a single energetic (or hard) gluon. The 
corresponding final state is expected to consist of three separated bundles of particles 
referred to as jets and originating from the quarks and gluon produced in the e"'"e~ 
interaction. These observables will be referred to as 3-jet observables. Some observables 
are defined such that they are only sensitive to 4-jet final states, i.e. those involving 
one gluon emission and one of the fundamental QCD processes listed above, and will be 
referred to 4-jet observables. 



2.2.1. Fixed order predictions The basic prediction in NLO for the normalised 
differential distribution of a 3-jet observable y is given by a power series in ds(Q) = 
«s(Q)/(27r): 

^ ^^ds(g) + ( (27r/3o ln(xD - 2) ^ + ^1 al{Q) . (16) 



cTtot dy <\y \^ ^ / dy J 

The energy scale Q at which the strong coupling is evaluated is generally identified 
with the physical hard scale of the process such as the cms energy of the e"'"e~ 
annihilation. The coefficient functions dA/dy and dB/dy correspond to the coefficients 
Ri and i?2 of equation (^Q) and their renormalisation scale dependence is given by 
equation © after replacing as by as and absorbing (27r)" in the coefficient i?„. The 
normalisation to the total hadronic cross section atot has been accounted for by the 
relation Utot = cro(l + 2as) [22]. The coefficient functions dA/dy and dB/dy may be 
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derived numerically for any suitable observable [23,24]. The derivation is done by Monte 
Carlo integration of the NLO QCD matrix elements for the production of up to four 
partons [22] over contours in phase space given by the definition of the observable y. In 
some cases the LO terms dA/dy have been obtained analytically, see e.g. [5]. 

A third order term ~ dC/dydg can be added to equation (fT?)|l once a numerical 
integration of the corresponding next-to-next-to-leading order (NNLO) QCD matrix 
elements is available to generate the NNLO coefficient function dC/dy [25,26]. The 
renormalisation scale dependence will be expressed as in equation (jH} while the 
normalisation to the total hadronic cross section must now consider cxtot = o"o(l + 
2q;s + 5.6368d;s) [27]. The complete result for the NNLO term is 

4vr/3o^ + (2-)^/5i^) + (2vr/3o Inx^ ^ ^ (17) 




dy ' ' ' 'dy J A* ' V ' - I'J dy 

-2 (^:^ + 2vr/?olnxJ^)- 1.6368^) ai(g) . 

For 4-jet observables NLO fixed order predictions are also available [28,29]. Their 
form is derived from equations (fT^ and (fTTj) by setting dA/dy = 0: 

2.2.2. NLLA predictions The NLO predictions described above work well in regions 
of phase space where radiation of a single hard gluon dominates. For configurations 
with soft or collinear gluon radiation from the initial quark-antiquark pair the Sterman- 
Weinberg criteria may not be fulfilled anymore and consequently new divergent terms 
appear. For observables y defined to vanish when no gluon radiation occurred, 
i.e. in the 2-jet limit, the typical leading behaviour of the cumulative cross section 
R{y)=jyi/ada/dy'dy' is: 

i?(y) ~ In^" ^ = a^L^" , (19) 

for each order n of the expansion in as with ln(l/?/) = L [30]. For y — > the simple 
NLO prediction will be unreliable as a^L"^ ^ 1 will not hold anymore. A systematic 
analytic resummation of the leading logarithmic terms ~ L^" and the next-to-leading 
logarithmic terms ~ L^""^ to all orders in as has been performed for a special class of 
observables [30]. Such observables exponentiate, which means that \n.R{y) ^ Lgi{asL) 
and Qi has a power series expansion in as-Z^ [30]. For the cumulative cross section the 
following representation is possible: 

R{y) =C{as)ny,as)+D{y,as) (20) 

oo 

C(as) = 1 + E Cn&^ 

n=l 
oo n+1 



lnE(|/,as)= Y.Y.Gnm&^L"^ 

n=l m=l 

= Lgi{asL) + 5'2(«s^) + • • • 
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where \im.y^QD{y,a^) = 0. The function gi contains the resummation of all leading 
terms ~ agL"+^ (LL) while g2 resums all next-to-leading terms ~ a^L"' (NLL). The 
description of equation (pn|l is expected to be valid in the region asL < 1 which goes 
further into the 2-jet region at small values of y than the NLO prediction bounded by 
asL"^ < 1. 

A general numerical method for performing the resummation for a large class of 
event shape observables has been presented in [31]. This makes resummation possible 
in cases where a purely analytical calculation is impossible or too difficult. 

2.2.3. Matched NLL A and fixed order predictions The NLO and NLLA predictions 
described above can be combined to yield a prediction which is valid in the 2- and 3-jet 
regions. In order to avoid double counting of terms ~ as and ~ a\ present in both 
NLO and NLLA predictions one has to identify and remove such terms from the NLLA 
prediction. Table El presents a comparison of the fixed order with the NLLA prediction 
for the quantity In Rijj) . 

Table 2. Comparison of NLLA and fixed order calculations. The dependence of the 
coefficient functions A, B, and C on y has been suppressed. 



LL 


NLL 


subleading 


non-log 


fixed order 








as 


= Aas 




G22alL^ 






^ {B^\A^)al 




Gss&lL^ 


Gaaai^" + G^ialL 




= {G-AB + \A^)al 


Lgi{asL) 


92{asL) 


+ ... 


+ ... 





The first two columns represent the resummation of leading and next-to-leading 
logarithmic terms while the equalities in the last column refer to summing each 
corresponding row. The coefficient functions A{y), B{y) and C{y) correspond to the 
cumulative cross sections, e.g. A{y) = dA/dy'dy' and analogously for B{y) and 
C{y). The normalisation of the coefficient functions is such that A^y^aax) = B{y^i^^) = 
C (^max) = where ymax is the maximum kinematically possible value of the observable. 
This normalisation generates all terms needed to account for the normalisation of the 
prediction to the total hadronic cross section atot- The Gnm coefficients in the first 
two rows are determined by expanding the functions gi and g2 in terms of as while 
subleading logarithmic terms can be determined numerically [30,32,33]. 

The resulting matched prediction is for the cumulative cross section: 

lni?(y) = A(y)ds + {B{y) - \A\y))al (21) 

+ Lgi{asL) + g2{a^L) 

— {G^L'^ + G'iiL)ds — {G^2L'^ + G22L)a\ . 

Renormalisation scale dependence and normalisation to the total hadronic cross section 
can be inserted into equation (pT|l by making the replacements Biy) — > B{y) + 
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2A(?/)(7r/3o Inx^ — 1) and G22 G22 + G'i227r/9o In as well as modifying §2 according 
to [30], equation (8). This form of matching is referred to as In (i?) -matching. Other 
forms of matching are possible with small differences in the treatment of higher order 
terms [30, 34, 35]. An important difference between the ln(i?)-matching and other 
matching schemes is that in the other schemes the coefficients Ci and G2 and in some 
cases G21 must be known explicitly. 

In order to force the NLLA terms to vanish at the upper kinematic limit |/max of 
an observable y the so-called modified matching is used [30,35]. The modified matching 
consists of replacing L = ln(l/y) by L' = ln(l + 1/y — l/?/max)- 

A possible NNLO term in the fixed order prediction as discussed above can be 
taken into account in the ln(i?)-matching scheme by adding the following term to the 
right-hand-side (RHS) of equation (j2T| : 

(C{y) - A{y)B{y) + ^A'{y) - G,,L' - G,,L'^ al . (22) 

Renormalisation scale dependence and normalisation to the total hadronic cross section 
can be considered by replacing Giy) = dG/dy'dy' by the correspondingly integrated 
equation (fTTj) and by the replacement G33 G'33 + G'2347r/5o In x^. For matching schemes 
other than the ln(i?)-matching the coefficients Gi,i = 1,2,3, G21 and in some cases 
Gsi, 2 = 2,1 must be available. 



2.3. Inclusive ohservahles 

In the definition of inclusive observables no requirements on the structure of the hadronic 
events are made, instead the cross sections or branching ratios for the production of 
hadronic final states are considered. The objects of interest are QCD induced corrections 
to electro-weak processes with hadronic final states. 

The total hadronic cross section in e"'"e~ annihilation ctq including 7 — Z° interference 
but without QCD corrections can be written as follows (see e.g. [5]): 

Ms) = ^[Qj- 2QfV,VfXi{s) + {Al + V^){A) + Vf)x2{s)\ , (23) 

where a is the QED coupling, Q/ is the fermion charge and the Aj and Vf are the 
axial and vector couplings of fermions to the Z°. These are given hy Af = Tf^^ and 
Vf = Tf^3 — 2Qfsm^Qw with T/^3 = +1/2 for neutral leptons and up-type quarks 
and Tj3 = —1/2 for charged leptons and down- type quarks with the electro- weak 
mixing angle Qw The functions xi = i^s^s — m|o)/((s — ?ti|o)^ + T'^om'^o) and 
X2 = /^^s^/((s — m|o)^ + r|om|o), K = \/2GFm\o / [IQua), take into account Z° exchange 
and its interference with photon exchange, respectively. 

The quantity -Re+e- = cr(e"'"e~ — >• hadrons)/(T(e"'"e" — * /U^/U^) is introduced to study 
the QCD corrections with suppressed electroweak effects. At values of the cms energy 
^/s far below the Z° resonance we have i?e+e- = A^c Hq Qq while for the hadronic partial 
decay width of the Z° one gets Rz = AcEg(^g + Vg)/{Al + V^^) with the number of 
quark colours Nc = 3. 
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The QCD induced corrections to -Rg+e- can be factorised such that one gets for 
< mzo [27,36,37]: 

i?e+e- ^ Re+e- (1 + 2ds + 5.6368d^ - 102.44d|) . (24) 

For the case of hadronic decays of r leptons produced in e"'"e~ annihilation the 
quantity Rr = r(r^ hadrons)/r(r^ z/T-z/ge") is considered. The QCD corrections 
lead to the following expression for Rr [27]: 

Rr = 3(1 + 2ds + 20.809d^ + 210.9d|) . (25) 

Section 16.21 gives more details about the QCD description of hadronic decays. 

2.4- Monte Carlo models 

Important tools for the analysis of hadronic e^e~ annihilation events are event 
generation programs based on the Monte Carlo method [38]. These programs simulate 
the production and development of individual hadronic events according to probability 
densities derived from the theory for the parton shower. In addition models are used to 
describe the transition from the parton to the hadron state, the so-called hadronisation. 
Brief reviews of the various programs can be found in [39, 40] while reviews of the 
underlying theoretical methods are e.g. [5,41]. Figure Ogives an overview of the various 
steps needed to simulate the production of an hadronic event. 

QCD calculations in the leading logarithmic approximation (LLA) form the basis 
of the parton shower simulation. The probability Pbc in the LLA for a parton a = q, g to 
branch to two partons be = qg, gg, qq depends on the strong coupling as{Q) evaluated 
at the energy scale Q of the branching and the kinematics of the branching process. 
The differential DGLAP equation [42-44] specifies how Pbc depends on the energy scale 
of the branching through the running of the strong coupling as'- 
HP r 

^ = CpasiQ') J pMdz (26) 

with t = 2lia{Q/ Aqcd) and z = Ef,/Ea, 1 — z = E^/Ea where Ei is the energy of parton 
i. Since the branching is understood to be a part of a tree-like Feynman graph the 
partons a, h and c will be virtual, i.e. off their mass shell. The splitting functions Pbc{z) 
take into account the different kinematics of the three possible branching processes 
q — > qg, g gg and g —>■ o^in the LLA. 

The task of a Monte Carlo parton shower algorithm is to first choose for a given 
parton at which value of the so-called evolution variable t the branching should take 
place and then to select a value of the energy fraction z. The branching probability Pbc 
given by equation (j2SI) is used to write the probability Pn.e.(tmax,t) that no emission 
takes place as the product of constant no-emission probabilities valid in small t intervals 
At. These intervals range from the maximum possible t value t^ax given by the process 
which produced parton a to the value t of interest: 

\nPn.e.{tmax,t) = ln(l - At) ^ -^^^ At (27) 
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dt 



in the limit At — > 0, — > cxd. Equation (|27p can be used to pick a value of t from a 
random number i? distributed evenly between and 1 with R = Pn.e.(tmax,t)- When the 
chosen value of t falls below a pre-defined value tmin corresponding to a cut on Q called 
Qo the parton cannot branch and is put on its mass-shell. This step will eventually stop 
the parton shower evolution when all remaining partons cannot branch anymore. 

Once a value of t is known the kinematics of the branching is determined by choosing 
values of z using the splitting functions Pbc{z). The azimuthal angle of the branching 
must be chosen as well to complete the configuration [5]. The exact definition of the 
energy scale Q of the branching has not been given yet, since the existing parton shower 
algorithms use different approaches, e.g. = or = [41]. 

The Monte Carlo parton shower simulation based on LLA QCD calculations shown 
so far does not include interference between soft but acolinear gluons produced in 
different branchings (colour coherence) [5,39,41]. It turns out that soft gluon interference 
can be treated in the LLA including subleading corrections by the introduction of 
angular ordering where the opening angles of successive branchings are required to 
decrease in addition to the values of t which decrease by construction. The destructive 
soft gluon interference is the QCD analogue to the Chudakov effect [45], where e+e" 
pairs produced at high energy are observed to generate less ionisation in the region 
where the two tracks of the e"'"e~ pair are close together. The intermediate photons 
involved in the ionisation processes cannot resolve the two individual charges of the 
nearby tracks thus leading to a suppressed ionisation rate. In the parton branching 
a ^ be a. subsequent soft but acolinear gluon cannot resolve the daughters be and thus 
it can be viewed as emitted by parton a. This can be implemented in the parton shower 
algorithm by the introduction of the angular ordering requirement. 

The parton shower algorithm, in particular with the angular ordering requirement, 
is not expected to describe effects of hard emissions well, since it is based on the LLA 
valid for soft and coUinear emissions. It is possible to correct the algorithms using 
the LO QCD matrix element for the hard emission process [46-49]. For theoretical 
consistency it is necessary to apply the correction to all branchings up to the hardest in 
the shower [47]. An alternative parton shower formulation based on the colour dipole 
approach [50] does not need such corrections since it uses the LO QCD matrix elements 
already at each branching. 

Further corrections using NLO QCD matrix elements to obtain a better description 
of multi-jet final states with hard and well separated jets are actively developed [51-55]. 

In the hadronisation phase hadrons are built from quarks and gluons left after 
the end of the parton shower by algorithms based on hadronisation models. In the 
hadronisation model it has to be specified how much of the energy and momentum of a 
parton is taken over by the hadron built from it. 

The so-called string model [48, 56] is based on the colour field line configuration 
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between a separating quark and antiquark. Naively, the colour field lines should look 
similar to electric field lines but the self interaction of gluons causes a contraction of the 
colour field lines into a narrow tube or string. As the quark and the antiquark move 
apart energy is stored in the stretching colour flux tube which can be expressed by a 
linear term in the potential VqcDir) ~ ^ + b-r between the quark and the antiquark. In 
this picture the colour flux tube or string is stretched until the potential energy suffices 
to create a new quark antiquark pair causing the string to break in two. This process is 
repeated as long as the partons have enough energy to support it. A gluon created in the 
parton shower will have strings attached to it following the colour flow and stretching 
to both the quark and the antiquark and will thus appear as a kink in the string. 

At the end of the string fragmentation quarks and antiquarks connected by strings 
arc treated as mesons. Production of baryons is included by the possibility to produce 
pairs of diquarks qq or qq when a string breaks or by the popcorn mechanism [57] 
where the baryons are made up from successively produced quarks or antiquarks. The 
detailed behaviour of the hadronisation is controlled by fragmentation functions and the 
parameter Qq specifying the parton energy at which the parton shower is terminated. 
A typical value for Qo is about 1 GeV safely above the Landau pole where perturbative 
QCD must fail. The fragmentation functions are probability densities f{z) for the 
fraction z of energy and longitudinal momentum {E+p{) the newly created hadron qiq2 
takes after a string connected to the quark qi broke to create a pair q2q2- Usually a 
fragmentation function 

/(^) ^ ^IzjLe^bmVz) (28) 

with ml = + pi and a and h free parameters is used. Alternatively a fragmentation 
function as proposed by Peterson et al. [58] 

_i__ii_)y' (29) 

may be used for the hadronisation of heavy quarks, i.e. b- and c-quarks. The parameter 
Cq is expected to scale like but has to be adjusted for each heavy quark flavour 

separately by comparison with experimental data. The momentum transverse to the 
string direction is given by a Gaussian distribution with the free parameter uq to adjust 
the width. 

The cluster hadronisation model [49,59-62] is based on the so-called preconfinement 
of colour in QCD [63,64], i.e. that partons in a shower build clusters of colour singlets 
with masses 0{Qq). The cluster model assumes that hadronisation is a local process 
and does not need fragmentation functions to describe hadronisation. In the model after 
the end of the parton shower gluons arc split into qq pairs and quarks close in phase 
space arc combined into colour singlet objects called clusters [60]. The masses of the 
clusters are 0{Qq) and are approximately independent of the hard scale Q. Clusters 
exceeding an upper limit on mass are further split into lighter clusters. The clusters are 
intermediate states in the hadronisation algorithm allowed to decay isotropically into 
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two hadrons, or only one hadron when the cluster is too light. Production of mesons 
and baryons is handled by creating either a quark pair qq or by creating a pair of 
diquarks qqqq. The cluster decay products are then randomly identified with hadrons 
fitting the quark contents with probabilities proportional to the spin degeneracy of the 
hadrons and the available phase space for the decay, i.e. the density of states. Only 
kinematically allowed decays are accepted and in the case of rejection the cluster decay 
algorithm starts again for the cluster in question until an allowed decay is found. 

After the hadronisation is finished the unstable hadrons are allowed to decay 
according to measured branching ratios and decay rates where known. Especially for 
hadronic decays of hadrons containing heavy quarks this information is not always 
available. In this case the heavy quark is allowed to decay weakly into quarks 
and the hadrons are created by the same hadronisation mechanisms as described 
above [48,49,65]. 

The next subsections summarise some of the most popular programs currently in 
use and briefly discuss the optimisation of the description of the data by the models. 

2.4- 1- PYTHIA The PYTHIA program [48] implements a parton shower algorithm 
essentially as described above with the invariant mass of the parton (virtuality) as the 
ordering parameter and angular ordering as additional constraint. The parton shower 
is combined with the string fragmentation model. 

2.4-2. HERWIG The HERWIG program [49] provides a parton shower similar to 
our description, but with the ordering variable replaced by the opening angle of the 
branching. This feature allows the inclusion of angular ordering in a clean way. 
Correlations between azimuthal angles of subsequent parton branchings due to gluon 
polarisation effects are also taken into account. The parton shower is terminated 
by setting minimal values for the invariant masses of quarks and the gluon. The 
hadronisation follows the cluster model including splitting of gluons into qq pairs after 
the parton shower stopped. In HERWIG clusters containing quarks produced in the 
parton shower decay such that the hadron carrying the quark flavour travels in the 
direction of the quark. Hadronisation of bb events can be controlled with additional 
parameters to take account of their special properties. 

2.4.3. ARIADNE The ARIADNE program [66] has a parton shower algorithm based 
on the colour dipole model [50] and with the relative momentum p\ of the branching as 
the ordering variable. The branchings are strictly ordered with 1 > 2 > P±,3 > ■ ■ ■ 
such that the angular ordering requirement is fulflUed. The parton shower terminates 
when all p\ are below a minimum value and the events are further processed by the 
string fragmentation model as in PYTHIA. 

2. 4-4- CO JETS The parton shower in the COJETS Monte Carlo program [67,68] is 
formulated using the LLA in a way so that a shower for each of the initial quarks or 
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antiquarks can be evolved separately. The first branching is constrained to the 0{a^) 
QCD matrix element for hard gluon radiation. There is no imposition of angular ordering 
in the parton shower, thereby neglecting the soft gluon interference effects. A limit on 
the energy of partons generated in the shower is used to stop the parton shower. The 
running strong coupling as is implemented using the one-loop expression. 

The hadronisation algorithm in CO JETS is based on the Field- Feynman model [69]. 
In this model each parton left after the parton shower has stopped is fragmented 
independently using an iterative algorithm. A quark qi is paired with an antiquark 
q2 from a newly created qq pair to build a hadron. The remaining quark q2 is then 
fragmented by the same procedure. After the fragmentation is finished some extra 
particles are added to ensure flavour conservation. The sharing of energy and momentum 
is specified by a fragmentation function with the Field- Feynman parametrisation f{z) = 
1 — a + 3a{l — z)"^ where a is a free parameter and z is the ratio of energy and momentum 
from the hadron and the quark qi. The models for the parton shower and hadronisation 
are more simple than those of the other programs; the COJETS program is now often 
used to study the effects of using more simple models. 

2.4-5. Tuning of Monte Carlo Models The free parameters of the models can be varied 
in order to achieve an optimal description of the data. The important free parameters 
to influence the behaviour of the parton shower are the value of the strong coupling 
constant as and the parameter to control termination of the shower. 

The string hadronisation is controlled by the parameters of the fragmentation 
functions for light or heavy quarks, see equations and (pUj) . and by the parameter 
setting the width of the Gaussian distribution of transverse momentum. Additional 
parameters regulate the production of baryons and strange hadrons. The cluster 
hadronisation model has free parameters to control splitting of heavy clusters into light 
clusters before these clusters are allowed to decay into hadrons. 

The LEP collaborations have studied the Monte Carlo models in detail and derived 
sets of parameters which optimise description of their data. The data consist of event 
shape observables, jet production rates, charged particle multiplicities and production 
rates of identified particles measured with the large event samples produced on the 
peak of the Z° resonance. The latest values of the model tuning parameters are given 
for ALEPH in [70], for DELPHI in [71], for L3 in [72], for OPAL in [73,74] and for SLD 
in [75]. 

2.5. Power corrections 

At energy scales Q ~ Aqcd perturbative evolution of as{Q) breaks down completely due 
to the Landau pole at Aqcd- It is therefore impossible to attempt pQCD calculations of 
soft processes at low values of Q without further assumptions. In the preceding sections 
Monte Carlo event generators were discussed which successfully model hadronisation 
effects. However, it is interesting to attempt to treat hadronisation in a more direct 
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way. 

In order to clarify which assumptions are needed a phenomenological model of 
hadronisation is considered. The longitudinal phase space model or tube model goes 
back to Feynman [76]. One considers a qq system produced e.g. in e"'"e~ annihilation in 
the cms system. The two primary partons move apart with velocity f/c ~ 1. In such 
a situation the production of soft gluons will be approximately independent of their 
(pseudo-) rapidity rj' = — log(tan(0j/2)), see figureEl 



A q 




Figure 3. Sketch of a qq system in rj' — kt space. 



The change to e.g. the observable Thrust 1 — T due to the production of a soft gluon 
at angle Qi with transverse momentum kt^i is A(l — T)j ~ fct,i/Qe~'^^'. This observation 
is generalised to write the soft or non-perturbative contribution to the mean value of 
the 1 — T distribution as follows: 

(l-T),p = /|$(A;0^-/e-I^W (30) 

The function $(fct) ~ tts(^t) is the distribution of soft particles in k^. The first integral 
in equation pO|l is summarised as ao/Q independent of the observable and the second as 
a constant ci_t dependent on the observable. The quantity ao ~ / «s(fct)dfct can only 
exist when as(fct) is identified with a non-perturbative strong coupling ag^(/i;t) which 
is assumed to be finite at low kt around and below the Landau pole. 

A more formal approach to the origin of soft contributions is the study of infrared 
renormalons, i.e. the divergence of asymptotic pQCD predictions due to the integration 
of low momenta in quark loops in gluon lines [77]. The infrared renormalon divergence of 
the C(ag) term in a pQCD prediction is factorial in n: r^a^^^ ~ (2/9o/p)""«!as'^^ With 
Stirlings formula to replace n\ one finds r^a^"*"^ ~ {2Poas/p)'^n"'e~^as- The convergence 
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of the series is optimal for n 



p/{2PQas). Based on this relation one finds 




QCD 



Q 



) 



(31) 



where the first order relation between as and Aqcd has been used. The result shows 
that infrared renormalon contributions to pQCD predictions scale like similar to 
the soft contribution studied in the tube model, see equation (I^Uj). 

The power correction model of Dokshitzer, Marchesini and Webber (DMW) 
extracts the structure of power correction terms from analysis of infrared renormalon 
contributions [78]. The model assumes that a non-perturbative strong coupling exists 
around and below the Landau pole and that the quantity ao(/ii) = l/fiiJo"^ a^^{kt)dkt 
can be defined. The value of /ij is chosen to be safely within the perturbative region, 
usually /ii = 2 GeV. A study of the branching ratio of hadronic to leptonic r lepton 
decays as a function of the invariant mass of the hadronic final state supports the 
assumption that the physical strong coupling is finite and thus integrable at very low 
energy scales [79]. 

The main result for the effects of power corrections on distributions F{y) of the 
event shape observables 1 — T, Mh and C is that the perturbative prediction FpT^y) is 
shifted [80-82]: 



where Cy is an observable dependent constant and P ~ M/ii/Q(ao(/^i) —C(s) is universal, 
i.e. independent of the observable [81]. The factor P contains the 1/Q scaling and the so- 
called Milan-factor M which takes two-loop effects into account. The non-perturbative 
parameter is explicitly matched with the perturbative strong couphng as. For the 
event shape observables Bt and Bw the predictions are more involved and the shape of 
the pQCD prediction is modified in addition to the shift [82]. For mean values of 1 — T, 
Mh and C the prediction is: 



For (Bt) and (Bw) the predictions are also more involved due to the modification of 
the shape of the distributions. 

3. QCD in e"'"e" annihilation 

Studies of QCD in e"'"e~ annihilation have been carried out for more than 30 years now. 
The earliest experiments were located at small colliders operating in the y/s range of 
about 1 — 2 GeV, ADONE, ACO and VEPP-2. A comprehensive review of the early 
experiments and results may be found in [83]. 

3.1. Accelerators and experiments 



F{y) = FMy - CyP) 



(32) 



(y) = {y)pT + 



(33) 



The development of storage rings with colliding beams of electrons and positrons was 
the foundation for the wealth of results we have today not only about QCD but also 
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about many other aspects of the Standard Model of particle physics. Table El collects 
the e"'"e~ storage rings and their experiments. The centre-of-mass energies available to 
the experiments in the laboratory system were increased by 2 orders of magnitude. 



Table 3. c+c colliders and experiments. 



Facility 



Location 



s [GeV] 



Experiments 



AGO [84] LAL Orsay ^ 1 

ADONE [87] INFN Frascati 1 - 3 

VEPP-2 [93] Novosibirsk 1 - L5 

CEA [94] Cambridge, MA 4 

SPEAR [96] SLAG Stanford 2-8 

PEP [101] SLAG Stanford 29 

DORIS [107, 108] DESY Hamburg 3-11 



CESR [115] 



Gornell, Ithaka 10-11 



PETRA [120] DESY Hamburg 12 - 47 



TRISTAN [125] KEK Tsukuba 50 - 64 



SLG [129] 
LEP [131] 



SLAG Stanford ^ 91 
GERN Geneva 88 - 209 



M3N [85,86] 
Boson [88], /xtt [89], 77 [90], 
772 [91], MEA [92] 
VEPP-2 
BOLD 
SLAG-LBL [97,98], 
MARK I [99], MARK II [100] 
MARK II [102], HRS [103], 
TPG/27 [104,105], MAG [106] 
PLUTO [109], DASP [110,111], 
LENA [112], DH(HM) [113,114] 
GLEO [116,117], 
GUSB [118,119] 
GELLO [121], JADE [122], 
MARK J [123], PLUTO [109], 

TASSO [110,124] 
TOPAZ [126], VENUS [127], 
AMY [128] 
MARK II [102], SLD [130] 
ALEPH [132,133], 
DELPHI [134,135], 
L3 [136], OPAL [137] 



The integrated luminosities of the corresponding data samples range from 
e'(l-100)/nb of the early experiments to 0(10 - 100)/pb at PETRA, PEP and 
TRISTAN. The experiments at the only recently decommissioned colliders SLG and 
LEP collected data samples corresponding to several hundred 1/pb on the Z° peak and, 
in the case of LEP, above the Z° peak up to cms energies of 209 GeV. 

The sizes of the data samples vary from C(IOO) events to 0(10 000) events per cms 
energy in experiments running at the colliders with cms energies below or above the TP 
peak. On the TP peak the LEP experiments accumulated about 5 ■ 10^ hadronic events 
each while the SLG experiment SLD collected (9(100 000) hadronic T^ decays. 
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3.2. Highlights of QCD before the LEP aera 

The first experimental studies of liadron production in e"'"e~ anniliilation wliere done 
in tlie early 1970s when evidence in support of the quark-parton model [138, 139] had 
come from deep-inelastic scattering experiments. 

A simple prediction of the quark-parton model was that the cross section for 
e"'"e~ hadrons should be large compared to the cross section for e"'"e~ — > fJ^^fi', 
because of the additional colour degree of freedom, see section Early measurements 
of -R = cr(e"'"e~ — >• hadrons)/(T(e+e~ — >• fJ'~^f^~) at y/s ~ 1 — 3 GeV indeed observed 
i? ~ 2 as predicted by the quark-parton model for u, d and s quarks [99] . However, the 
measurements were difficult to interpret quantitatively due to the presence of resonances. 

The quark-parton model also predicts that hadron production in e"'"e~ annihilation 
at sufficiently high energy should show a pattern of two jets of hadrons recoiling against 
each other. The direction of the hadron jets is given by the direction of the produced 
quarks and should thus follow the expected distribution of the angle 9 between the quark 
and the beam direction for the production of two particles with spin 1/2 [140-142]: 

-^-l + cos^e . (34) 
dcos H 

The first evidence for jet structure was reported 1975 in [143] using data recorded at 
y/s = 6.2 and 7.4 GeV by the SLAC-LBL magnetic detector at the e+e~ storage ring 
SPEAR. Hadronic events were analysed with the event shape observable Sphericity [141]. 
Distributions of Sphericity values and angular distributions of jet directions given by the 
Sphericity axis were compared with jet model and phase space Monte Carlo simulations. 
This evidence was further supported in [144-146]. 

QCD as the field theory of quark interactions predicts the existence of gluons as 
intermediate gauge bosons analogous to the photons of QED. It was crucial for the 
acceptance of QCD that effects directly connected with gluon activity could be observed 
in hadronic e^e~ annihilation events. 

The first indirect evidence of the existence of gluons came from so-called direct 
decays of the Y(1S) resonance into hadrons. In [147] the PLUTO collaboration using 
data recorded on and off the Y(1S) at the DORIS storage ring at DESY studied event 
shape observables like the Thrust (see section In direct decays of the Y(1S) the 

Thrust distribution was found to agree with Monte Carlo simulations based on Y(1S) 
decays mediated by three gluons while simple phase space models were ruled out. A 
study of the angular distribution of the Thrust axis w.r.t. the beam direction provided 
evidence that the gluons are vector particles, i.e. have spin 1. 

The first direct evidence for the existence of gluons was provided by the observation 
of planar hadronic events with a clear 3-jet structure by the experiments at the PETRA 
e"'"e~ collider at DESY [148-151]. The cms energies of the collisions ranged from 17 to 
32 GeV. In these studies events were classified as deviating from a 2-jet configuration 
based on the event shape observables Sphericity, Thrust or Oblateness. Examination 
of the energy fiow in the event plane [148-150] showed a pattern consistent with LO 
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QCD based Monte Carlo models including gluon radiation. Distributions of event shape 
observables Oblateness, Thrust or Planarity were measured in [149,151] and were found 
to be reproduced by Monte Carlo models including gluon radiation at leading order. 
In [151-153] the angular and energy distribution of the lowest energy jet are found to 
agree with the expected behaviour for radiation of energetic vector gluons. 

The first measurements of the strong couphng constant as were based on Monte 
Carlo simulations based on LO QCD and data from the PETRA experiments recorded at 
^/s 30 GeV [151, 152, 154-156]. The simulations generated qqg final states according 
to LO QCD matrix elements and employed a simple Field-Feynman fragmentation 
model. The value of as together with parameters to adjust the fragmentation model was 
varied until the description of event shape observable distributions or 3-jet production 
rates was optimised. The results were as(30GeV) = 0.19 ± 0.04 corresponding to 
as(mzo) = 0.15;ni 

Improved determinations of as using NLO QCD calculations turned out to be 
difficult to achieve [122]. For event shape or jet observables the early NLO calculations 
gave differing results and the dependence on the fragmentation models resulted in large 
errors. For inclusive quantities like the R-ratio Re+e- the sensitivity was limited by 
uncertainties of the available data. Global averages shortly before the start of LEP 
were as(mzo) = 0.11 ±0.01 [157] based on most available measurements and as(mzo) = 
0.112 ± 0.015 from event shape and jet observables alone [158]. These determinations 
had a relative uncertainty of about 10% where the dominating uncertainties came from 
the use of fragmentation models and missing higher orders in the perturbative QCD 
calculations. 

The non-abelian nature of QCD becomes manifest in its property of asymptotic 
freedom and in the existence of the triple gluon vertex (TGV) . The first evidence for 4- 
jet structure was found by studying event shape observables constructed to be sensitive 
to non-planar multi-jet events [159]. The first experimental evidence for the running 
of the strong coupling as, i.e. for asymptotic freedom, was based on comparing data 
from PETRA runs at several cms energies between 22 and 47 GcV in terms of 3-jet 
fractions determined with the JADE recombination jet algorithm [160, 161]. The 3-jet 
fractions were shown to correspond closely to the underlying parton structure assumed 
in QCD with unimportant dependence on fragmentation models and thus to reflect the 
running strong coupling. In [162] the evidence was confirmed and in addition angular 
distributions measured with 4-jet events were found to prefer a Monte Carlo model 
based on 0{al) QCD including the TGV over an abelian model. 

Dedicated studies of properties of hadronic e"'"e~ events showed that fragmentation 
models using independent fragmentation of the hard partons could not describe the 
data completely, in particular the production of soft particles between jets, the so-called 
string-effect, see e.g. [122,163,164]. This reduced the set of models to those which used 
the Lund string model or the cluster hadronisation model [165]. It was shown that the 
string effect can be explained in QCD by coherent emission of soft gluons [166]. In 
addition it was found that production rates of multi-jet events were better reproduced 



CONTENTS 



26 



by models using a parton shower algorithm [167] based on QCD in the LLA, with 
further improvements when the parton shower is matched to the 0(as) QCD matrix 
element [46, 161]. 

To summarise, the quark-parton model had been established by the earliest e^e~ 
annihilation experiments. Fundamental predictions of QCD, the field theory of strong 
quark interactions, such as the existence of gluons and their properties, had been 
successfully tested. However, the limited precision of calculations and the lack of 
fundamental understanding of the fragmentation process did not allow experimental 
tests of the theory with uncertainties better than about 10% as reflected in the errors of 
the combined values of as(^zo) shown above. In particular the final experimental proof 
of asymptotic freedom had not yet been demonstrated although there was already good 
positive evidence. Many detailed studies of the properties of the hadronic final states in 
e"'"e~ annihilation had shown that fragmentation models implementing a parton shower 
with coherence effects together with the Lund string or the cluster hadronisation model 
could describe the data most successfully. 

4. Jets and event shapes 

One of the most prominent and important features of hadronic final states produced in 
e"'"e~ annihilation is the jet structure, i.e. the presence of a small number of collimated 
groups of particles recoiling against each other, see also section El Quantifying the 
structure of hadronic final states allows for direct comparison of experimental data with 
predictions by QCD and thus for detailed and stringent tests of the theory. Many 
attempts have been made to quantify the jet structure; these schemes fall in two 
categories: jet clustering algorithms and event shape observables. The jet algorithms 
and prescriptions to calculate event shape observables must be infrared and collinear 
safe in order for perturbative QCD predictions to be possible, see section [21]. 

The following sections discuss commonly used jet clustering algorithms and event 
shape observable definitions. Experimental results will be presented and compared to 
model predictions. Direct tests of asymptotic freedom in QCD and measurements of as 
using jet or event shape observables will be discussed in separate sections. 

4-1. Jet observables 

In jet clustering algorithms one tries to group the particles of the hadronic final state 
such that the jet structure, which is often clearly visible, is replicated. A popular 
algorithm was introduced by the JADE collaboration [160,167], often referred to as the 
JADE algorithm. 

The first ingredient of the clustering algorithm is the definition of the distance yij 
in phase space between any two particles i,j in the hadronic final state. The distance 
measure is defined as 




(35) 
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where Mij is the invariant mass between two particles i,j and E'vis = Hi Ei is the total 
visible energy of all particles in the hadronic final state. The second ingredient is the 
prescription for the combination of two particles i,j into a cluster or jet fc; in the case 
of the JADE algorithm this is done by adding the 4-momenta p: pk = Pi +Pj- The 
algorithm proceeds recursively by combining the pair of particles with the smallest yij, 
removing the particles i,j from and adding the jet k to the event until yij > Ucnt for all 
remaining y^. The remaining jets are the result of the clustering algorithm. 

Many variations of this scheme exist differing in the calculation of yij and the 
combination prescription, see e.g. [168]. The original JADE algorithm (JADE EO) 
defines Mf^ = 2EiEj{l — cos%) taking all particles as massless. The Durham 
algorithm [169] defines y^j = 2mm{Ef , E — cos 9ij)/ E^^^, i.e. the distance is given by 
the relative transverse momentum between one of the particles w.r.t. the other. 

The JADE EO and Durham algorithms were shown to be preferable, because i) 
they have little dependence on hadronisation models compared to other algorithms and 
ii) they were found to have comparatively small theoretical uncertainties [168, 170]. 
The Durham and JADE EO algorithms have hadronisation corrections of about 5 % at 
^/s = mzo [168]. The Cambridge algorithm [171] is a variant of the Durham algorithm 
where some soft particles or jets are excluded from further jet clustering in order to 
reduce hadronisation effects. However, in [172] it was shown that hadronisation effects 
for 2-, 3- and 4-jet fractions are actually larger for the Cambridge than for the Durham 
algorithm in large regions of ycnt, but that these effects cancel for the more inclusive 
observable given by the average multiplicity of jets per event. 

The events are usually analysed by studying the fraction of e.g. 3- or 4-jet events 
as a function of the jet resolution parameter ycut- Alternatively, the values of ?/cut at 
which the number of jets in the events change e.g. from 2 to 3 are used to compute 
differential distributions [173]. In this case the observable is analogous to the event 
shape observables discussed below. 

A large number of results on jet production in e"'"e~ annihilation has been published; 
reviews of earlier work can be found in [174-176]. Recent and comprehensive studies of 
jet production may be found in [72, 177-179]. 

Figure m shows the fractions of 2-, 3-, 4- and 5-jet events using the Durham (left) 
and JADE (right) algorithms measured by L3 using the highest energy LEP 2 data at 
an average ^/s = 206 GeV [72]. The data are compared with predictions by Monte 
Carlo models and reasonable agreement is found. 

Figure El displays in (b) the 3- and in (c) the 4-jet fraction determined with the 
Durham algorithm at ^/s = 35, 91 and 189 GeV [177]. The 35 GeV data stem from 
a reanalysis of data from the JADE experiment while the 91 and 189 GeV data are 
from OPAL. These figures show that the same jet fractions are found for smaller values 
of ?/cut at increasing ^/s. In the regions of large j/cut w.r.t. the peaks this effect is a 
consequence of the running of as- The small decrease in the maximum 3-jet and 4-jet 
fractions with increasing ^/s is also due to the running of as. The data are compared 
to various Monte Carlo models tuned to OPAL data at i/i ~ mzo. Good agreement is 
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Figure 4. Jet production rates with the Durham (left) and JADE (right) algorithms 
measured at various values of t/cut by L3 [72]. The 2-, 3-, 4- and 5-jet production 
rates corrected for experimental effects are shown as points with error bars; the lines 
indicate the corresponding Monte Carlo model predictions as indicated on the figures. 



found for all models except COJETS which fails to describe the data at high energy. 




Figure 5. 3- and 4-jct production rates with the Durham algorithm measured by 
JADE at Vs = 35 GeV (triangles) and OPAL at = 91 GeV (points) and 189 GcV 
(squares) [177]. The data are corrected for experimental effects. The predictions by 
Monte Carlo models are indicated by the hues for PYTHIA (soUd), HERWIG (dashed), 
ARIADNE (dotted) and COJETS (dash-dotted). 

While the jet production rates are generally well described by the usual Monte Carlo 
event generators it has been observed that the detailed kinematics of 4-jet events are 
not always successfully modelled [180-182]. The kinematics of 4-jet events are studied 
using angular correlations between the jets, e.g. the Bengtsson-Zerwas angle [183] 
Xbz = ^{\pi X P2], [Ps X P4]), where the pj, i — 1 . . .4 are the momentum vectors of 
the four jets after energy ordering. Comparing LEP 1 data with simulations showed 
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deviations of up to about 20% while the data can be well described by LO {0{a\)) or 
NLO {0{a\)) calculations corrected for hadronisation effects [180,181,184,185]. 

The jet structure of hadronic events in e"'"e~ annihilation has also been studied using 
a cone based algorithm as commonly used at hadron colliders [177,186]. An advantage 
of this study is that the properties of the jets may be compared more directly to results 
from hadron collider experiments. 

4-2. Event shape ohservahles 

The event shape observables avoid direct association of particles to jets and calculate 
instead a single number which classifies the event according to its jet topology. Generally 
the observables are constructed such that a value of zero corresponds to an ideal 2- 
jet event consisting of only two back-to-back particles, a small value corresponds to 
a realistic 2-jet event while increasingly larger values indicate the presence of one or 
more additional jets. A large number of event shape observables have been proposed 
and studied. We will concentrate here on observables which are infrared and collinear 
safe and which have the most complete QCD predictions including contributions from 
resummed NLLA calculations. 

Thrust T, Thrust major and minor Tmaj. and Tmin.: These observables are de- 
fined by the expression [187, 188] 



where Pi is the momentum of particle i in an event. The thrust axis nx is the vector 
n which maximizes the expression in parentheses. A plane through the origin and 
perpendicular to ht divides the event into two hemispheres Hi and H2- A value of 
T = 1 corresponds to an ideal 2-jet event; therefore the thrust observable is often 
used in the form 1 — T. Based on the thrust axis ht the thrust major Tmaj. and the 
thrust major axis ut^^^ are defined by equation ()36|) with the constraint nx-n = Q. 
The thrust minor Tmin. and the thrust minor axis ut^^^ are defined by the expression 
in parenthesis of equation with the constraint [ut x ut^^^ ] = nr^^^ ■ 
Heavy and Light Jet Mass Mh and M^: The hemisphere invariant masses are 
calculated using the particles in the two hemispheres Hi and H2. We define 
Mh [189, 190] as the heavier mass and Ml as the lighter mass, both normalised 

to y/s. 

C- and D-parameter C and D: The linearised momentum tensor G"'^ is defined by 




(36) 




The three eigenvalues Xj of this tensor define C [22] through 



C = 3(AiA2 + A2A3 + A3Ai) 
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where Oij is the angle between particles i and j. The D-parameter is defined by 



Jet Broadening observables Bt, Bw and Bn'- These are defined by computing 
the quantity 



for each of the two event hemispheres, if^, defined above. The three 
observables [191] are defined by 



where Bt is the total, By/ is the wide and B^^ is the narrow jet broadening. 
Transition value between 2 and 3 jets y^^: The Durham jet algorithm (see sec- 
tion I4.1|l is applied to cluster the particles of an hadronic event into two or three 
jets. The value of Ucnt at which for an event the transition between a 2-jet and a 
3-jet assignment occurs is called [173]. 

We will refer to event shape observables in general by the symbol y. The observables 
Tmin., Ml, D and B^ belong to the class of 4-jet observables, because for these we have 
in the cms y > only for events with at least four partons in the final statef. We also 
note that some of the observables take the whole event into account while others depend 
only on one selected hemisphere of the event; the observables Mh-, Ml, Bw, B^ and 
1/^ fall into the latter category. 

Figure El presents measurements by OPAL at average y/s = 91, 133, 177 and 
197 GeV of the event shape observables (left) and D (right) representing 3- and 
4-jet observables [182]. The data are corrected for experimental effects and compared to 
predictions by Monte Carlo models. One finds that the data for yf^ are well described 
by the models at all energies. The data for D are also well described by the models at 
high energies ^/s > mzo where the statistical and experimental errors are comparatively 
large. However, at ^/s = 91 GeV and large values of D significant deviations between 
models and data occur; for HERWIG there are also discrepancies at ~ 10^^. This 
observation could be related to the unsuccessful modelling of the kinematics of 4-jet 
events discussed above. 

Another way of studying the structure of hadronic events in e"'"e~ annihilation is 
to calculate moments of the event shape observable distributions: 



fDue to momentum conservation in the cms these observables have j/ = for back-to-back two 
parton as well as for planar three parton configurations. 



D = 27A1A2A3 




Bt = Bi + B2, Bw = max(5i, B2) 



and Bn = min(i?i, B2) 
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Figure 6. Distributions of event shape observables (left) and D (right) measured 
by OPAL [182]. The data are corrected for experimental effects. The high energy 
LEP 2 data collected at several energy points are combined into two samples at average 
-/s = 177 and 197 GeV. The data sets are separated vertically by multiplicative factors 
as indicated for clarity. The lines show the predictions of Monte Carlo models as 
indicated on the figures. The lower parts of the figures present the deviations between 
data and Monte Carlo model predictions at ^/s = 91 and 197 GeV. 

where |/max is the maximal kinematically allowed value of the observable. The moments 
always sample all of the available phase space with a different weighting depending on 
the value of n; moments with small values of n are dominated by 2- and 3-jet events while 
for larger values of n multi-jet events become more important. Figure [7| (left) shows 
the first five moments of the observables 1 — T, Mh and C measured by OPAL [182] 
at y/s = 91 to 197 GeV. The data are compared to predictions by Monte Carlo models 
and good agreement is found, except for HERWIG which shows deviations from the 
moments of 1 — T and C increasing with n. 

The reanalysis of data from the JADE experiment made it possible to compare data 
for event shape observables at y/s < nizo with the predictions of current Monte Carlo 
models [192,193]. Figure [7| (right) presents as an example measurements of distributions 
of 1 — T at ^/s = 14, 22, 35 and 44 GeV. The data from the JADE experiment are 
corrected for experimental effects and the presence of e"'"e~ —>■ bb events corresponding 
to about 9% of the event samples. The predictions of Monte Carlo models tuned by 
OPAL to LEP 1 data are superimposed and found to agree well with the data; the 
exception is COJETS which deviates from the data in particular at low i/s. 
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Figure 7. (left) Moments n = 1, . . . , 5 of event shape observables 1 — T, Mh and C 
measured by OPAL [182]. The data are corrected for experimental effects. The lines 
show the predictions by Monte Carlo models as indicated. The lower parts of the figure 
present the deviations between data and Monte Carlo models at ^/s = 91 and 197 GeV. 
(right) Distributions of 1— T measured using JADE data and corrected for experimental 
effects at y^=14 to 44 GeV [192,193]. The lines or crosses show predictions by Monte 
Carlo models as indicated. The set of figures on the RHS present the correction factors 
Ki (lines) and their uncertainties (shaded bands), and the differences 5i between data 
and Monte Carlo model predictions together with the errors of the data. 

^.5. Experimental tests of asymptotic freedom 

Jet rates and event shape observables are well suited to investigate experimentally one 
of the most important predictions of QCD, namely asymptotic freedom of the strong 
coupling, see sections |21 and El The large range in cms energies ^/s provided by the LEP 
data and also the earlier data from previous e"'"e~ colliders gives direct access to effects 
caused by the running strong coupling constant. 

Following earlier studies [160, 161] the L3 collaboration measured the fraction of 
3-jet events using the JADE EO algorithm at ?/cut = 0.08 [72]. This measurement 
is especially suited to study asymptotic freedom, since the hadronisation corrections 
as predicted by Monte Carlo models are known to be small (less than 10%) and fairly 
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model independent [168,170]. The L3 data and results from other experiments spanning 
the range ^/s = 22 to 206 GeV are shown in figure |H1 (left) and compared with the 0{al) 
perturbative QCD expectation with asifrizo) = 0.12. The prediction clearly agrees with 
the data. 




Figure 8. 3-jet fraction /a with the JADE EO algorithm at ycut — 0.08 (left) and 
1/(1 — T) (right) measured as a function of the cms energy [72, 194]. The data are 
compared with predictions by QCD shown by the solid lines. 



Figure |H1 (right) presents the results of a similar analysis by DELPHI [194] using 
the mean values (1 — T) of the 1 — T distributions measured at energies ranging from 
below 10 GeV to 202 GeV. The figure shows the quantity 1/(1 — T) vs. ^/s on a 
logarithmic scale, because from equations ((Zj) and (fTHjl a logarithmic dependence in LO is 
predicted. This prediction is well confirmed by the data. The solid line shows the QCD 
prediction in NLO in a special formulation independent of a particular renormalisation 
scheme [14,195] (see section ll.5.2|) which is also seen to suppress non-perturbative effects. 



4-. 4- Measurements of as 

4.4-1 ■ 3-jet observables Early measurements of as before the start of LEP using event 
shape or jet observables were briefiy summarised in section El The first analyses by 
the LEP collaborations used 0(ag) QCD calculations and already a set of Monte Carlo 
event generators including QCD coherence effects and tuned to the precise LEP data, 
see e.g. [175] for a review. An average value of as(mzo) = 0.120±0.006 was determined, 
where the error is dominated by theoretical uncertainties estimated by variation of the 
renormalisation scale in the QCD calculations. Compared to the earlier determinations 
of as(^zo) the error of the combined LEP measurements was reduced by a factor of 
more than two. This improvement was mostly due to the use of consistent hadronisation 
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models as well as reliable 0(ag) QCD calculations [23]. The error was also seen to be 
consistent with the scatter of results from individual observables; this is an important 
cross-check of the methods for estimating the errors. A reanalysis of low energy data 
at = 35 GeV from the JADE experiment using the methods developed for LEP 
obtained as(mzo) = 0.122l|5:g{]^ consistent with the LEP result [196]. 

A significant improvement in the determinations of as based on event shape or jet 
observables was the introduction of resummed calculations matched to the fixed order 
predictions (0(q;s)-|-NLLA) discussed in section |21 These calculations allow the data to 
be described over larger regions of phase space in the 2-jet region and have generally a 
weaker dependence on the renormalisation scale [34]. A combination of values of aa,{mzo) 
obtained from analyses at -y/i ~ mzo using (9(ag)+NLLA calculations was given in [197] 
as as(mzo) = 0.121 ± 0.006. Thus the total error was not improved compared to 0{a'^) 
based analyses but the result is considered to be more reliable theoretically due to the 
reduced dependence on the renormalisation scale. 

Recently the LEP experiments have published their final results on determinations 
of as based on data at i/i = mzo up to the highest LEP 2 energies [72, 178, 182, 198]. 
These analyses use a consistent set of observables: 1 — T, Mh, Bt, Bw and C and in 
addition by ALEPH and OPAL. All analyses employ the most recent calculations 
for the resummed NLLA calculations [30-33] and the modified ln(i?)-matching scheme 
and thus the results may be compared directly. 

Figure IHl (left) shows as a representative example the results of the analysis by 
ALEPH using the observable 5^ at a/s = 91.2 to 206 GeV [178]. The data are corrected 
for experimental effects and the curves present the 0(ag)+NLLA QCD prediction using 
the modified In (i?) -matching corrected for hadronisation. The fits have been performed 
over ranges indicated by the solid lines while the dashed lines are an extrapolation 
of the theory using the fitted results. The agreement between the theory including 
hadronisation corrections and the data is good in the fitted regions and becomes worse 
in the extreme 2-jet and multijet regions where the theory is not expected to be able to 
completely describe the data. 

In order to study asymptotic freedom the individual results on as from the most 
recent LEP analyses [72, 178, 182, 198] are combined using a common method, see 
e.g. [178,182,199,200]. The values of as to be combined are evolved to the common scale 
of the combination if necessary. Then the average value as is computed by minimising 

= {a-as)-V-' -{a-asf (37) 

w.r.t. to as- The individual measurements as^j are collected in the vector a and V is 
the covariance matrix of the measurements. The final expression for the average value 
is 

"s = X!^i«s,i , Wi = ^^^^ . (38) 

The correlations between experimental, hadronisation and theoretical uncertainties are 
generally not well known, in particular between different experiments, and thus models 
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Figure 9. (left) Distributions of Bt measured by ALEPH at several energy 
points as indicated [178]. Superimposed are fits of the QCD prediction corrected for 
hadronisation effects. The fits have been performed in the ranges shown by the solid 
lines while the dashed lines show an extrapolation of the fitted predictions, (right) 
Combined values of as based on the analyses by the LEP experiments or using JADE 
data (see text) are shown. The total (inner) errors bars show total (experimental) 
uncertainties. The lines and shaded band show the evolution of as based on the value 
of as(mzo) given on the figure. 

for the correlations have to be constructed. In these models three levels of correlation 
Pij between the errors (Xj and (Tj of two measurements i and j may be assumed: i) no 
correlations pij = 0, ii) partial correlation pij = min((jj, cXj)^/ {(XiCrj) or iii) full correlation 
pij = 1. The treatment of the different systematic uncertainties is as follows: 

Experimental uncertainties These are added to the diagonal and off-diagonal 
elements of V using one of the three possible models for their correlation. 

Hadronisation uncertainties These are added only to the diagonal elements of V. 
The hadronisation uncertainty is evaluated by repeating the combination with 
individual as^i values obtained with different hadronisation models. 

Theoretical uncertainties These are also added only to the diagonal elements of 
V. The theoretical uncertainty is evaluated by varying all individual values as,i 
simultaneously within their theoretical uncertainties and repeating the combination. 

This procedure suppresses the presence of negative weights Wi which otherwise appear 
when large off-diagonal terms are included in V [200]. Such negative weights are 
mathematically correct but difficult to interpret as a physical result. We therefore 
choose a treatment of the correlations of the generally large hadronisation and theoretical 
uncertainties which avoids negative weights as in [178,182]. 

In a first step the individual results for each observable and experiment are 
combined in the energy ranges 130 < ^/s < 136 GeV, 161 < y/s < 183 GeV and 
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189 < -y/i < 207 GeV corresponding to average centre-of-mass energies y/s = 133, 177 
and 197 GeV. In these combinations the statistical errors are taken as uncorrelated 
and the experimental errors as partially correlated. OPAL publishes results directly at 
^/s = 177 and 197 GeV and also at points from smaller energy ranges. Our averages 
reproduce the OPAL results within 1% and with consistent uncertainties. The single 
DELPHI result from ^/s = 183 GeV is simply evolved to 177 GeV. 

The sets of results for each observable and experiment are further combined 
within each experiment. The statistical correlations between observables as obtained 
by ALEPH and OPAL at y/s = 206 GeV are taken from [200]. The experimental 
uncertainties are again taken as partially correlated. Our values for as are shown in 
table m The corresponding combined results from OPAL have been calculated using the 
same method except for a more detailed treatment of hadronisation uncertainties. Our 
combined results again reproduce the OPAL results within 1% with consistent errors. 

Table 4. Results for average values of as from LEP experiments [72, 178, 182, 198]. 
The columns show the results for as{\/s), the statistical (stat.), experimental (exp.), 
hadronisation (had.) and theoretical (theo.) uncertainties. The last two rows give 
average values of as{mza) of all LEP 2 data from y/s = 133 to 197 GeV and of all 



LEP data. 

Q;s(m2o) ± stat. ± exp. ± had. ± thco. 

ALEPH 0.1201 ± 0.0002 ± 0.0007 ± 0.0020 ± 0.0050 

DELPHI 0.1216 ±0.0002 ±0.0021 ±0.0024 ±0.0056 

L3 0.1211 ±0.0008 ±0.0014 ±0.0035 ±0.0050 

OPAL 0.1193 ±0.0002 ±0.0007 ±0.0024 ±0.0047 

LEP 0.1204 ±0.0002 ±0.0006 ±0.0025 ±0.0051 

as (133 GeV)± stat. ± exp. ± had. ± theo. 

ALEPH 0.1156 ± 0.0028 ± 0.0008 ± 0.0014 ± 0.0048 

L3 0.1133 ±0.0023 ±0.0018 ±0.0024 ±0.0048 

OPAL 0.1091 ±0.0033 ±0.0033 ±0.0016 ±0.0036 

LEP 0.1124 ±0.0017 ±0.0011 ±0.0018 ±0.0043 

as (177 GeV)± stat. ± exp. ± had. ± thco. 

ALEPH 0.1106 ± 0.0024 ± 0.0008 ± 0.0010 ± 0.0040 

DELPHI 0.1093 ± 0.0035 ± 0.0020 ± 0.0017 ± 0.0043 

L3 0.1077 ±0.0020 ±0.0014 ±0.0018 ±0.0045 

OPAL 0.1074 ±0.0010 ±0.0011 ±0.0009 ±0.0032 

LEP 0.1081 ± 0.0013 ± 0.0009 ± 0.0013 ± 0.0040 

as(197 GeV)± stat. ± exp. ± had. ± theo. 

ALEPH 0.1075 ± 0.0012 ± 0.0009 ± 0.0008 ± 0.0036 

DELPHI 0.1083 ± 0.0013 ± 0.0021 ± 0.0017 ± 0.0042 

L3 0.1124 ±0.0007 ±0.0013 ±0.0014 ±0.0046 

OPAL 0.1074 ±0.0010 ±0.0011 ±0.0009 ±0.0032 

LEP 0.1085 ± 0.0006 ± 0.0007 ± 0.0011 ± 0.0038 

as('7i2o) ± stat. ± exp. ± had. ± theo. 

LEP 2 0.1200 ± 0.0007 ± 0.0010 ± 0.0016 ± 0.0048 

LEP all 0.1201 ±0.0005±0.0008±0.0019±0.0049 



Applying the same procedure to the reanalysed data from the JADE experiment 
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collected at ^/s = 14 to 44 GeV [201] yields the results shown in table EH". The 
analysis of the JADE data was performed using the same event shape observables, 
QCD calculations and essentially the same Monte Carlo programs to calculate the 
hadronisation corrections as for the LEP analyses and therefore the results may be 
compared directly. 

Table 5. Results for average values of as from JADE. The last row gives the average 
value of asimzo) from all JADE data. 





Q!s(Vs) i stat. ± exp. ± had. ± theo. 


14 


0.1704 ± 0.0044 ± 0.0041 ± 0.0148 ± 0.0117 


22 


0.1514 ± 0.0040 ± 0.0030 ± 0.0106 ± 0.0084 


35 


0.1430 ± 0.0010 ± 0.0022 ± 0.0076 ± 0.0079 


38 


0.1388 ± 0.0035 ± 0.0030 ± 0.0056 ± 0.0070 


44 


0.1299 ± 0.0020 ± 0.0031 ± 0.0057 ± 0.0057 


as('7izo)± stat. ± exp. ± had. ± thco. 


JADE 


0.1203 ±0.0007±0.0017±0.0053±0.0050 



The average values of ag determined at each ^/s as shown in tables ^ and El are 
presented in figure IHl (right). The data are compared with the QCD prediction for the 
running of as calculated in NNLO based on the average value as{mzo) = 0.1211±0.0021 
(see section [7j) which has been determined using only NNLO QCD calculations. The 
agreement with the data is excellent and provides strong evidence for asymptotic 
freedom. 

Measurements of the first five moments of the six event shape observables discussed 
above and as shown partially in figure [3 (left) have been used by OPAL to determine the 
strong coupling constant using only NLO (0(0^)) QCD predictions [182]. The NLLA 
predictions break down at very small values of the observables y <^ 1 and also become 
unreliable close to the kinematic limits ymux and thus cannot be used to derive moments. 
The divergences of the NLO predictions for y —>■ are not a problem, since these are 
integrable and stable predictions for the moments are obtained. After demanding that 
the ratio of the NLO to the LO terms in the perturbative QCD predictions be less 
than 0.5 results from 17 observables are considered: (1 — T), {C), (Bt), {By^) and 
(1/23) with n = 1, . . . , 5, and (M|}) with n = 2, ... ,5. The final combined result is 
as(mzo) = 0.1223 ± 0.0005(stat.) ± 0.0014(exp.) ± 0.0016(had.)+|;;|!J5J(theo.) in good 
agreement with the results from distributions discussed above. This analysis constitutes 
a valuable cross check of the determination of as from 3-jet observables, since with 
moments the complete phase space is used. 

4.4-2. 4-j^i observables The study of 4-jet observables for measurements of as(mzo) is 
promising, because the sensitivity of these observables to the strong coupling is doubled 
compared to 3-jet observables. The relative change in as induced by a change of a 3-jet 

fThe statistical errors at = 14 and 22 GeV were estimated using the statistical error at 38 GeV. 
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observable O3 ~ as is given by Aas/o^s = AO3/O3 while for a 4-jet observable O4 ~ 
we have Aag/ag = A04/(204). 

The ALEPH collaboration has performed a measurement of as{mzo) using the 
fraction of 4-jet events R4 determined using the Durham algorithm in hadronic 
decays [181]. The QCD prediction is known in 0{a\), i.e. LO is 0{a\) and NLO 
radiative corrections are C(ag) [28,202]. These calculations are matched with existing 
NLLA calculations. 

The measurement of Oisirnio) based on these predictions is presented in 
figure ^1 [181]. The uncorrected data for R4 measured with 2.5 ■ 10^ hadronic TP 
decays using the Durham algorithm are shown as a function of the jet resolution 
parameter ?/cut- The line shows the fit of the (9(ag)+NLLA (R- matching) theory 
corrected for hadronisation and detector effects. The fit is performed using only 
the bins within the fit range as indicated where the total correction factors deviate 
from unity by less than 10%. The result of the fit for fixed renormalisation scale 
parameter = 1 is as(mzo) = 0.1170 ± 0.0001(stat.) ± 0.0022(syst.) using the 
conservatively estimated systematic uncertainty from [181]. The systematic error has 
contributions from experimental (±0.0009), hadronisation (±0.0010) and theoretical 
(±0.0017) uncertainties and the fit is seen to describe the precise data fairly well. 

A similar analysis by DELPHI [203] uses only the 0{a\) (NLO) calculations 
with experimentally optimised renormalisation scale and finds as('^zo) = 0.1175 ± 
0.0005(stat.)±0.0010(exp.)±0.0027(had.)±0.0007(theo.) with = 0.042 for the closely 
related Cambridge algorithm [171] and as{mzo) = 0.1178 ± 0.0036 with = 0.015 
for the Durham algorithm. A study from OPAL using R4 measured with LEP 1 and 
LEP 2 data found as(mzo) = 0.1182 ± 0.0003(stat.) ± 0.0015(exp.) ± 0.0011(had.) ± 
0.0018(theo.) [204]. Measurements of R4 with JADE data at y/s between 14 and 
44 GeV found as a preliminary resuh as(mzo) = 0.1169 ± 0.0004(stat.) ± 0.0012(exp.) ± 
0.0021(had.) ± 0.0007(theo.) [205]. The various measurements are in good agreement 
with each other. All measurements using R4 discussed in this section have comparatively 
small uncertainties and belong to the group of the most precise determinations of as- 

4-5. Alternative approaches to soft and hard QCD 

4-5.1. Tests of power corrections Distributions and mean values of event shape 
observables measured in e"'"e~ annihilation at many cms energy points between around 
14 and more than 200 GeV have been used recently by several groups to study 
power corrections [72,178,194,198,206]. In all of these studies pQCD predictions in 
(9(Q;g)±NLLA are used for event shape distributions. For the mean values only NLO 
(C(ag)) predictions are employed as explained in section 11.4. II The pQCD predictions 
are added together in both cases with the power correction predictions. The resulting 
expressions are functions of the strong coupling ag(mzo) and of the non-perturbative 
parameter ao. The complete predictions are compared with data for event shape 
distributions or mean values corrected for experimental effects. For the mean value 
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Figure 10. The 4-jet fraction R4 determined using the Durham algorithm is shown 
as a function of j/cut without corrections. Superimposed is a fit of ^^(aD+NLLA QCD 
(sohd hne) corrected for hadronisation and detector effects [181]. 



of the 4-jet observable D-parameter power correction predictions have been tested as 
well [72]. 

Figure im (left) shows the results of fitting the data for 1 — T distributions measured 
at cms energies from 14 to 189 GeV. The data from experiments at y/s < rrizo are 
corrected for the effects of e"'"e~ —>■ bb events [206]. The data in the fitted regions are 
well reproduced by the theory (solid lines). The dotted lines present extrapolations 
outside of the fitted regions using the fit results for as(mzo) and ao(2 GeV) which 
describe the data reasonably well. The agreement between theory and experiment is 
similar for other observables. 

Figure ITTl fright) presents as solid lines the results of fits to measurements of mean 
values of various event shape observables at y/s = 12 to 206 GeV [178]. The data are well 
described by the fitted theory. The dotted lines show the results of fits where the power 
corrections have been replaced by corrections derived from Monte Carlo simulations. 
The Monte Carlo corrections lead to slighty different theoretical predictions. 

Table IHl collects averaged results for as(mzo) and ao(2 GeV) from the various 
analyses [72,178,194,206]. The averages are as given by the authors and are based 
upon differing sets of observables, definitions of systematic uncertainties and averaging 
methods. It is therefore not useful to compute global averages from these results. There 
is good agreement between the various analyses of mean values or distributions of event 
shape observables. 

Figures ^1 (left) and (right) give an overview over results from the analysis of 
distributions [178] (left) and mean values [72] (right). The results are shown as points 
with one-standard-deviation error ellipses. The results for as('^zo) of fits to distributions 
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Figure 11. (left) Distributions of 1 — T. The vertical axis has arbitrary units and 
the data points have been scaled for clarity. The solid lines show the result of the fit of 
the combined pQCD and power correction prediction. The dotted lines represent an 
extrapolation of the result [206]. (right) Cms energy dependence of the mean values of 
several event shape observables. The solid lines show the result of the fit of combined 
0{a'^) QCD calculations with power corrections while the dashed lines indicate results 
obtained with Monte Carlo based corrections [178]. 



Table 6. Rc'sults of po\v(>r correction analysis using t he DMW model. 

ALEPH [178] DELPHI [194] L3 [72] [2061 





as(mzo) 


0.111 


0.108 




0.111 


distri- 




±0.005 


±0.002 




±0.004 


butions 


ao(2 GeV) 


0.50 


0.55 




0.58 






±0.10 


±0.03 




±0.09 




as(mzo) 




0.121 


0.113 


0.119 


mean 






±0.006 


±0.006 


±0.003 


values 


ao(2 GeV) 




0.47 


0.48 


0.49 








±0.08 


±0.06 


±0.06 



arc found to be smaller compared to the results of fits to the same data using Monte 
Carlo based hadronisation corrections; as{mzo) ~ 0.111 compared to as{mzo) ~ 0.120. 
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as(Mz) 



Figure 12. Summaries of results for as(wzo) and ao(2 GeV) for distributions 
(left) [178] and mean values (right) [72]. The error ellipses show one standard deviation 
uncertainties (39% CL). The shaded band on the left plot gives the combined result 
for Q!g(mzo) using Monte Carlo based corrections from the same analysis. The bands 
on the right plot give the unweighted averages of as(m-zo) a-nd ao(2 GeV) shown in 
this figure. 

There are systematic differences between the DMW power correction model and the 
hadronisation corrections derived from Monte Carlo models leading to differences in 
the values of as(mzo) of Aas(mzo) ^ 0.009 which are significantly larger than the 
hadronisation uncertainties quoted in section ll.4.1l as derived from the published results. 
In the analysis by ALEPH [178] a direct comparison is made using the same data and 
a consistent average relative difference Aas("^zo)/tts("^zo) — 10% is found. The results 
for ao(2 GeV) from p = Mfj, and Bt are fairly large compared to the averages. For 
the observable Mh this effect stems from neglecting particle masses [76,194]. 

Figure Uni shows a summary of results for asirriio) and q;o(2 GeV) [72,178,194,206]. 
One observes reasonable agreement between the individual measurements. The non- 
perturbative parameter ao appears universal within about 20% as expected [81] while 
the values for as{mio) are compatible with world averages, e.g. [207] a;s(mzo) = 
0.1184 ±0.0027. 

4 -5. 2. Studies of renormalisation schemes The ideas about choosing optimal 
renormalisation schemes (RS) discussed in sect ion 1^ . 1 . 1 1 have been tested experimentally. 
In [197] the predictions for from the three methods for the hadronic width of the TP 
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Figure 13. The figure shows a summary of results for q;s(?7t-zo) and q;o(2 GeV) from 
power correction studies [72,178,194,206]. The resuhs labelled PPTSL are from [206]. 



Rz = r(Z° hadrons)/r(Z° — > leptons) were obtained for NLO QCD. Since for Rz the 
QCD prediction is known in NNLO the resulting values of as(mzo) from scale optimized 
NLO QCD may be compared with the corresponding result from NNLO QCD. It is found 
that for this observable the three methods agree well, but the resulting value of as(mzo) 
is significantly lower by Aas{mzo) — 0.003 compared to the values determined with the 
NNLO calculations. 

In [208] the DELPHI collaboration measured many event shape observables using 
precise LEP data from the Z° peak and compared the data to NLO QCD predictions; a 
similar analysis is discussed in [209]. In fits where the renormalisation scale parameter 
Xfj, and thus effectively the RS was allowed to vary in addition to as{mzo) consistent 
results within the errors for as(^z") were found resulting in an unweighted average 
of Q!s(^zo) = 0.1170 ± 0.0025 where the error is given by the variance of the set 
of measurements. Within the fitted ranges of the distributions no evidence for a 
dependence of on the observable was found, implying that a single value of is indeed 
sufficient. As a comparison the corresponding fits with = 1 yielded an unweighted 
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average and variance of aQ,{mio) = 0.123±0.015; i.e. not using experimentally optimised 
scaled results in a significantly increased variance. Similar observations were made 
in [210]. However, these results are not supported by the findings of [209] where for a 
similar set of observables as in [208,210] as{mzo) = 0.127 ± 0.008 for fits with = 1 
and as(mzo) = 0.117 ± 0.007 for fits with free was observed. 

Applying the scale or RS setting methods to fits of limited ranges of distributions 
is in principle problematic, because more than one energy scale is involved. However, 
since it is known experimentally that a single value of is sufficient to describe the 
data the effects of the additional energy scales introduced by the observable values may 
be assumed to be small. In [208] the ECH and PMS methods generally lead to similar 
results for as{mzo), consistent with the theoretical observation that their predictions for 
the optimal scale are closely related, see e.g. [15]. The resulting values of as(^zo) are 
also consistent with each other. The weighted averages of the set of observables with the 
PMS or ECH methods are as(^zo) = 0.115 ± 0.004 with xVd-o.f. ~ 1. For the BLM 
method x^/d.o.f. = 29/13 is found indicating less consistency among the results. The 
observation that values of as{mzo) obtained from theoretically optimised RS tend to be 
smaller than purely experimental values may be explained by noticing that in the PMS 
method (which is close to the ECH method) a stationary point and usually a minimum 
of as{mzo,Xfj) is chosen as the optimum. 

By extending the ECH method it turns out to be possible to relate QCD predictions 
for different physical observables to each other without explicit dependence on the 
RS [14,17,211]. The QCD prediction for an inclusive observable R{Q) depending on 
only one energy scale Q is given by the generalised /3-function: 



with b = 27r/9o- The pi are RS independent but can be derived from the corresponding 
RS dependent coefficients. Equation may be solved analogously to the /3-function 
of as in terms of a process dependent scale parameter Apt- Thanks to the Celmaster- 
Gonsalves relation [212] such parameters for different observables may be related to each 
other or e.g. to the scale parameter of the running strong coupling as{p) in the MS 
renormalisation scheme: 



The variables are Ci = 2tt(3i/ (3q and ri = Bji/{2A^ where and are the LO and 
NLO coefficients for the observable calculated in the MS scheme. This procedure may 
be referred to as renormalisation group improved perturbative QCD (RGI QCD) [10]. 
Non-perturbative effects scaling like 1/Q for most observables may be included with an 
additional free parameter Kq [213]. 

The DELPHI collaboration has studied mean values of several event shape 
observables and in addition distributions of the event shape observables EEC and JCEF 
using the RGI QCD method as explained above [194]. The analysis also uses data from 
lower energy experiments to cover values of ^/s from 22 GeV to 200 GeV. The results for 



Q— = -bR\l + piR + p2R^ + ...) = bp{R) 



(39) 




(40) 
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Q;g(mzo) from 12 observables are consistent with each other within their uncertainties, 
as shown in figure (left). The variance of the individual values of a^irriio) is about 
Aas(mzo) ~ 0.002 while the variance of as(^zo) results of corresponding fits based on 
standard NLO QCD combined with power corrections is observed to be about twice 
as large. The values of the power correction parameter Kq are found to be small and 
consistent with zero within their uncertainties. Figure El (right) shows as examples the 
results of RGI QCD fits without power corrections, i.e. Kq = 0, for the mean values of 
Tmaj., 1 - T and C. 




0.09 0.095 0.1 0.105 0.11 0.115 0.12 0.125 0.13 0.135 0.14 VsrGeVl 



Figure 14. (left) The figure summarises results for asimza) from RGI QCD fits 
with power corrections. The error bars include statistical and experimental systematic 
uncertainties. The solid lines and grey bands indicate the unweighted mean and 
variance based on mean values of event shape observables. (right) The figure shows 
fits of RGI QCD without non-perturbative corrections to mean values of Tmaj., I — T 
and C as functions of (full lines). As a comparison the dashed lines present the 
corresponding prediction in the MS RS [194]. 



The results of the RGI QCD fits are interpreted to imply that a dominating part 
of the power corrections or hadronisation corrections in standard NLO QCD fits in the 
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MS RS can be absorbed in the RGI QCD. In [194] a systematic shift in as(?7izo) of 
about 2% between fits of RGI QCD with and without power corrections was observed 
and this shift is conjectured to be an estimate of the true size of missing higher order 
contributions or power corrections. 

In [214] the RGI QCD method is applied to extract a value for as('^zo) from the 
branching ratio of hadronic to leptonic decays of r leptons. The result turns out to be 
in good agreement with the current world average. See section 16.31 for a discussion of 
the extraction of as from hadronic r decays. 

4-6. Running b quark mass 

The QCD analyses discussed so far have been made assuming massless partons in the 
QCD predictions. The effects of massive quarks on the perturbative QCD predictions are 
known in 0(a;g) as well [215-217] (and also in NLLA [218]). The mass of the heavy quark 
mb identified with the b quark at LEP energies becomes an additional parameter of the 
theory and is subject to renormalisation analogously to the strong coupling constant. 
As a result the heavy quark mass mb will depend on the scale of the hard process in 
which the heavy quark participates. In leading order in the MS renormalisation scheme 
the running b-quark mass is mb((5) = Mb(as(Q)/7r)^^/^^, where Mb is the so-called pole 
mass defined by the pole of the renormalised heavy quark propagator (see section !^. 1.3|1 . 
The running mb is known to four- loop accuracy [19]. Using low energy measurements 
of mb(mb) one finds that mb(mzo) ~ 3 GeV. 

The mass of the b quark mb is experimentally accessible in hadronic Z° decays, 
because i) about 21% of the decays are Z° bb and these events can be identified 
with high efficiency and purity, ii) event samples of O(IO^) events are available 
and iii) observables like the 3-jet rate R3(2/cut) in the JADE or Durham algorithm 
can have enhanced mass effects going like m'^/m^o/ycut [219]. In the experimental 
analyses [220-224] the ratio B3 = R^s/R's, R''3(R^3) are the 3-jet rates in b quark 
(light quark) events, is studied which reduces the influence of common systematics. 

Figure HSl (left) shows Bad/cut) measured by OPAL [221] before and after detector 
and hadronisation corrections for the JADE EO and the Durham algorithm. One 
observes that the effect of the b quark mass is about 5% on B3 with reasonable 
uncertainties as indicated by the bands. The effect goes in opposite directions for 
JADE EO and the Durham algorithm, because in the JADE type algorithms based on 
invariant mass to cluster jets the large b quark mass causes an enhancement of B3 while 
in the Durham algorithm the reduced phase space for hard gluon emission dominates 
thus reducing B3 [222]. 

Figure El (right) gives an overview over presently available measurements of 
mb(mzo) [224]. These are compared with low energy measurements and the QCD 
prediction for the running of mb(Q). The predicted range of mb(mzo) is in good 
agreement with the measurements. The measurements of mb(mzo) of [221-224] are 
averaged using the same method as for averaging values of as (see section I4.4.1|) . 
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Figure 15. (left) Distributions of B3 determined using the Durham or JADE EO jet 
finding algorithms before and after all corrections. The solid lines in the lower row of 
figures give the uncertainties on B3. The vertical dashed lines indicate the values of 
?/cut where the QCD predictions are calculated [221]. (right) Measurements of TOb(Q) 
at low and high energies Q compared with the QCD expectation for the running of 
mb(Q) [224]. 



Table 7. Values of my,{mz(>) in imits of GeV from various experiments with 
symmetrised systematic uncertainties. The last row shows the average determined 

as described in the text. 

Experiment rrih (tozo )± stat. ±syst . ± had. ±theo. 

ALEPH [223] 3.27± 0.22 ±0.22 ±0.38 ± 0.16 

DELPHI [224] 2.96±0.18±0.14±0.19± 0.12 

OPAL [221] 2.67±0.03±0.22±0.25±0.19 

SLD [222] 2.56±0.27±0.33±0.84±0.44 

average 2.90±0.12±0.29±0.20±0.14 



The individual results are shown with symmetrised systematic uncertainties in table [3 
The systematic uncertainty of the OPAL result was split into an experimental and a 
hadronisation uncertainty according to the average ratio of these uncertainties obtained 
from the individual jet algorithms. We assume the statistical errors to be uncorrelated 
and the experimental and hadronisation uncertainties to be partially correlated. The 
theory uncertainties are evaluated by simultaneous variation of the input values within 
their theory uncertainties. The average is shown in the last column of table [7| 
Combining all errors we find mb(mzo) = (2.90 ± 0.31) GeV, which may be compared 
with the average of low energy measurements mb(mb) = (4.24 ±0.11) GeV [225]. The 
difference becomes mb(mzo) — mb(mb) = (1.34 ± 0.33) GeV which is non-zero by four 
standard deviations thus providing strong experimental evidence for the running heavy 
quark mass in QCD. 
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5. Particle production 

The direct observation of liadron production in e"'"e~ anniliilation offers unique ways to 
test predictions of QCD. The subject of multiplicities of hadrons is covered in detail 
in [226]. 

5.1. Low energy particles 

Figure El (left) shows the production rates for inclusive and for identified charged 
hadrons (tt^, K^., p) or hadrons decaying to final states with charged particles {K^, 
0, A) as a function the particles scaled momentum fraction x = 2^/1/5 where p is the 
particle momentum [227]. The momentum spectra of the various particles are fairly 
similar but their production rates are strongly suppressed with particle mass. The 
production rates for many particle species are summarised in [18]. 

The study of the observable ^ = ln(l/x), where x is the scaled momentum or energy 
fraction defined above, allows to focus on the production mechanisms of low-momentum 
particles corresponding to large values of ^. Figure Uni (right) shows the ^ distributions 
measured by OPAL using all LEP 1 and LEP 2 data together with data from previous 
experiments [228]. A similar analysis is shown in [230]. Figure IT7I (left) presents recent 
measurements by BES of ^ at = 2.2 to 4.8 GeV [231]. 

On both figures the typical hump-backed plateau [232] of the distributions is clearly 
visible. The distribution vanishes for small values of ^ corresponding to large x for 
kinematic reasons while the fall for large ^ (small x) is due to destructive interference 
of coherent gluon radiation. Essentially, low energy gluons cannot resolve individual 
coloured partons anymore and are thus suppressed since an e"'"e~ annihilation event 
is colourless. This explanation in the framework of perturbative QCD rests on the 
assumption that hadronisation is on average a local process. This assumption is referred 
to as local parton hadron duality (LPHD) [64,232]. Assuming LPHD predicted spectra 
calculated for partons may be directly compared with measured spectra derived from the 
observed hadrons where only the relative normalisation is determined by experiment. 
The calculation in the modified leading-log approximation (MLLA) [233-235] is shown 
as solid lines in figure IT7I (left) while the related NLLA QCD prediction [229] is shown 
as lines in figure EJ (right). 

The position of the peaks ^* of the ^ distributions is seen to depend on the cms 
energy. This effect is predicted by MLLA QCD [235] and shown as the line on figure IT7I 
(right). The fit was performed using only the low energy BES data and the result agrees 
well with the data up to LEP 1 energies. 

Measurements of ^ and ^* have also been made for identified hadrons (vr^, K^'", 
p. A) as well as for identified quark fiavours [130] and dependence on hadron mass and 
quark flavour has been observed. Distributions of C, at different y/s may be related 
directly via perturbative evolution equations in the MLLA [236]. 
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Figure 16. (left) Spectra of scaled momentum x for identified particles as indicated 
and for all charged particles measured by SLD on the Z° peak [227]. (right) Spectra of 
^ = log(l/a;) measured by OPAL using LEP 1 and LEP 2 data as well as lower energy 
data from TASSO [228]. The curves show fits of a NLLA QCD prediction [229] to the 
data. 

5.2. High momentum particles and scaling violation 

The fragmentation of partons produced in the hard interaction e+e~ qq(fl') into 
hadrons is investigated using the scaled momentum of hadrons x. In the quark parton 
model, i.e. without strong interactions, one expects the distributions of x measured 
at different cms energies to be identical. This prediction of a scaling behaviour is 
similar to the situation in deep inelastic scattering (DIS) of leptons on hadrons where 
XBj = Q^/{2m{E — E')) with transferred 4-momentum Q, parton mass m, and lepton 
energies E and E' before and after the scattering. Distributions of XBj are independent 
of in the quark parton model. 
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Figure 17. (left) Measurements of ^ = log(l/a;) spectra at 2.2 < ^/s < 4.8 GeV by 
BES [231]. (right) Measurements of the peak position of ^ distributions [231]. 



The scaling behaviour is modified or broken by strong interaction processes [237]. 
The probabihty for radiation of gluons from quarks depends on the energy scale of the 
quark and therefore due to the running of as fewer hard gluons are produced at high 
energy. The scaled momentum x spectra are expected to be softer at high energies. 
This is formally described by the theory with the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) evolution equations, see [238] and references therein. The normalised 
distribution of x is described as the differential cross section for production of a 
hadron with scaled momentum x: 

— -T—=f H Cf{z,as{nR),Xf,)Df{-,fi^)— . (41) 

The coefficient functions C/(x, Q;s(/iR), x^) correspond to the probability to obtain a 
parton / with scaled momentum x, where / is either one of the five active quark 
fiavours u, d, s, c, or b or a gluon g. The parameter x^ = fiji/Q, with Q = ^/s, 
gives the renormalisation scale /xr where the coefficient and fragmentation functions are 
evaluated. 

The fragmentation functions Df{x,fi) correspond to the probability to obtain a 
hadron with momentum fraction x from a parton / at the momentum scale of the hard 
interaction fi. The coefficient functions are known in NLO [240,241]. The fragmentation 
functions are non-perturbative and thus cannot be calculated; they are the e"'"e~ analog 
of parton density functions in DIS. However, the rate of change (evolution) of the 
Df{x,fi) with changing hard scale /x is predicted by the DGLAP equations: 

r = E / Pi/(^,«s(/iR),x^)A(-,/i)— (42) 
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Figure 18. (left) Ratio of tlie inclusive cliarge(i particle x spectra 1 / sigmahda^ / dx 
nieasureci by ALEPH at ^/s ~ 91 GeV to data from TASSO measured at y/s = 
22 GeV [239]. The lines present the result of a fit of the NLO QCD prediction (see 
text); the dashed shows the prediction with constant quark flavour composition, (right) 
Measurements of the normalised cross section for charged particle production in bins 
of a; as a function of ^/s [228]. The solid lines show the result of a fit of the NLO QCD 
prediction (see text) while the dashed show an extrapolation of the fit result. 



where the sphtting functions Pij are also known in NLO QCD, see [238] and references 
therein. With these predictions one measurement at a reference scale Qq may be used 
to predict the x distribution at another scale Q with the strong couphng constant as as 
a free parameter. 

Figure^! (left) shows the ratio of x distributions measured by ALEPH at ^/s = mzo 
and by TASSO at ^/s = 22 GeV [239]. The smaller (larger) fraction of hadrons at large 
(small) X at LEP 1 energies compared to ^/s = 22 GeV is clearly visible. The solid line 
indicates the result of a fit of the NLO QCD prediction for this observable to several 
measurements of x distributions at different energy scales. Figure UHl (right) presents 
the measurements by OPAL using LEP 1 and LEP 2 data together with measurements 
by previous experiments [228]. The solid lines show the QCD fit as explained above 
with as = 0.1182 fixed and variable parametrisations of the fragmentation functions. 
The description of the data by the fit is adequate within the uncertainties. The result 
of a fit with as(^«zo) as an additional free parameter is shown in table |H1 

Table IHl collects results for as (mzo) from analysis of scaling violation in x spectra 
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Table 8. Values of as{mio) determined from analyses of scaling violation of x 
spectra in e'^e^ annihilation. Experimental errors also include statistical uncertainties 
and uncertainties from parametrisation of non-perturbative fragmentation functions. 
The calculation of the average value is explained in the text. 
Experiment as(™zo)i expt. ± theo. 



ALEPH [239] 0.126 ± 0.007 ± 0.006 

DELPHI [242] 0.124 ± 0.007 ± 0.009 

OPAL [228] 0.113 ± 0.005 ± 0.007 

average 0.1192 ±0.0056±0.0070 



from e"'"e~ annihilation to hadrons. All analyses are based on the same NLO QCD 
calculations for the coefficient and splitting functions and thus the results may be 
compared directly. We observe good agreement between the three values of asijnio) 
within the experimental uncertainties. The average value of as{mzo) from analyses of 
scaling violation in e"'"e~ annihilation is derived assuming the experimental uncertainties 
to be partially correlated while the theoretical uncertainties are accounted for by varying 
the input values simultanously within their theoretical errors. The result is shown 
in the last row of table |H1 and is in good agreement with other determinations of 
asirriio) discussed in this report. The weights of the individual results in the average 
are 0.35 (ALEPH), 0.15 (DELPHI) and 0.50 (OPAL). This determination of as(mzo) 
does not have explicit hadronisation uncertainties, because the inclusive measurement 
of hadron production results in a supression of non-perturbative effects with l/Q"^ [243]. 
Remaining hadronisation uncertainties enter as uncertainties of the parametrisation of 
the fragmentation functions Df and are thus part of the experimental errors. 

Measurements of x distributions for identified hadrons [227] have also been studied 
in global QCD analyses using perturbative evolution [244, 245] and good agreement 
between data and theory was found. The scaling violations of identified light and 
heavy (b) quark and gluon jets using LEP 1 and LEP 2 data has been measured and 
compared to NLO QCD predictions in [246]. The QCD predictions for x spectra have 
been observed to be generally in good agreement with the data. 



5.3. Longitudinal and transverse cross section 

The measurement of the polar angle dependence of hadron production provides insight 
into the spin structure of the process e"'"e~ — > hadrons [238]. The normalised double 
differential cross section for hadron production in the process e"'"e~ ^ h + X is written 

as 

^ "^^""^ FT{x)l{l + cos^e)+FL{x)^sm^e + FA{x)lcose (43) 



o"tot da;d cos 9 8 4 4 

where a'* is the cross section for producing a hadron h, x is the scaled momentum fraction 
of the hadron and 6 is the angle between the hadron and electron beam directions. The 
fragmentation functions (EEs) are given by Fp{x) = l/atotdcp/da;, P = T,L,A. The 
first two terms stem from the polarisation states (transverse or longitudinal) of the 
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intermediate vector boson (7 or Z") while the third term comes from parity violation of 
the electroweak interaction. 

After integrating out the 9 dependence only the Ft{x) and Fl{x) terms remain. 
As a result of energy conservation the integral 



1 /-i 1 da^ 



X- 



dx 



{x){n) 



(44) 



2 Jo (Ttot dx 2 ' ■ ' ■ CTtot CTtot 

where (n) is the average hadron multiplicity and 1/2 / xFp{x)dx = ap, P = T,L. The 
transverse and longitudinal cross sections ctt and ctl fulfil ctt + ctl = crtot- 

Measurements of the FFs Fp{x), P = T,L,A for charged hadrons have been 
performed by ALEPH, DELPHI and OPAL using large samples of hadronic 
decays [239,247,248]. The associated cross section ratios (TT/crtot and crL/crtot for all 
hadrons are derived from the Fp for charged hadrons assuming that the ratio of charged 
to neutral hadrons is identical for ctt and (Jl- This assumption is supported by studies 
of simulated events [247,248]. A re-analysis of JADE data gave results for ctt and (Jl 
at yi = 36.6 GeV [249]. 

Figure Ca (left) shows the FFs determined by DELPHI and OPAL [248]. The 
analyses agree well with each other and a Monte Carlo model prediction (JETSET) 
is also found to agree with the data. The longitudinal FF F^ is suppressed w.r.t. to 
the transverse FF Ft reflecting the preferred polarisation of the intermediate gauge 
bosons. The presence of F^ is understood as due to the possible presence of gluons in 
the flnal state. The results for ctl are shown in figurelTTIl (right) [249] and compared with 
the JETSET Monte Carlo model prediction for partons and hadrons. The difference 
between the two lines for partons and hadrons indicates the presence of non-perturbative 
corrections. These power corrections have been calculated [243, 250] and were found to 
agree with the data [249]. 
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Figure 19. (left) Fragmentation functions Ft and Fl measured by DELPHI and 
OPAL [248]. The lines show the prediction of the JETSET Monte Carlo model, 
(right) Measurements of (Jh/crtot by DELPHI, OPAL and JADE compared with parton- 
and hadron-level predictions by the JETSET Monte Carlo model and a NLO QCD 
prediction [249] (see text for details). 
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The extraction of as(mzo) from the data can be done using the NLO prediction 
for the observable ai^/a^ox. [240,248]. In order to obtain the best possible measurement 
we average the results for ai^/atot from ALEPH, DELPHI and OPAL. The ALEPH 
measurements of F^^x) and Ft{x) [239] are used to extract a value for ai^/crtot after 
estimating the missing contributions at < x < 0.008 using the ALEPH measurement 
of (nch.) = 21.91±0.22 [251] and the DELPHI measurement oi Fl{x)/{Fl{x) + Ft{x)) = 
0.286±0.021 for < a; < 0.01 [248]. The estimate uses the relation /J FL{x)+FT{x)dx = 
{rich.)- The result is (crL/crtot)ALEPH = 0.0561±0.0003(stat.)±0.0023(syst.). The average 
of the ALEPH, DELPHI and OPAL results is computed assuming the experimental 
uncertainties to be uncorrelated with the result UL/crtot = 0.0558 ± 0.0021. The value 
of as(mzo) is extracted based on the NLO QCD prediction [240,248] combined with a 
hadronisation correction using JETSET derived from [247]. We find 

ag(mzo) = 0.1169 ± 0.0035(exp.) ± 0.0018(had.) ± 0.0072(theo.) . (45) 

The hadronisation uncertainty is given by the difference of values for Q;g(mzo) when 
HERWIG instead of JETSET is used to derive the hadronisation correction [247]. The 
theoretical uncertainty is evaluated by varying the renormalisation scale parameter in 
the range 0.5 < x^ < 2.0 according to equation ^ and is consistent with [248]. 

The connection of the asymmetric FF Fa{x) with parity violation in the electroweak 
interaction involving quarks has been investigated in [252]. A precise measurement of 
Fa{x) could be converted into a measurement of the electroweak mixing angle sin 6w 
assuming that hadrons at high x are likely to carry a quark produced in the e"'"e~ 
annihilation. This prediction has been found to agree with data in [247, 248] but the 
precision of the measurements was insufficient to extract a value for sinOw. 

5.4- Fragmentation of b quarks 

The fragmentation process with heavy quarks is of special interest, because the large 
mass of the heavy quark provides a natural cutoff in perturbative QCD calculations (see 
e.g. [253,254]). As a consequence reliable predictions of the perturbative part of the 
fragmentation process become possible. Assuming the b quark fragmentation process 
to be universal, precise measurements of the fragmentation of b quarks into b-flavoured 
hadrons with e"'"e~ annihilation data can have consequences for the interpretation of 
hadron collider data for B meson production [255]. 

Figure 1201 (left) presents as a representative example the measurement by SLD of 
the X spectrum of B hadrons produced in e"'"e~ annihilation at ^/s ~ mzo [256]. The 
variable x = xb = '^^Es/y/s corresponds to the scaled momentum variable discussed 
above, but due to the large B hadron masses the B hadron energy Eb is used. Other 
recent measurements are [257,258]. 

The measurement of the x spectrum of B hadrons is model dependent, because 
correcting the measured B hadron energies for experimental effects uses Monte Carlo 
simulations with models for the b quark fragmentation. The first step in Monte Carlo 
simulations is the parton shower which corresponds to a description of the perturbative 
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Figure 20. (left) Normalised x spectrum of B liadrons produced in e^e^ annihilation 
on the Z° peak [256]. (right) Difference of charged particle multiplicity between light 
and b quark events 5h\ measured by OPAL and other experiments at various cms 
energies [258]. The data are compared with several predictions as explained in the 
text. 



part of the transition of a heavy quark into a observable heavy hadron. After the 
parton shower has stopped hadronisation models such as the Lund string model describe 
the formation of hadrons from the remaining quarks and gluons of the parton shower. 
The heavy quark fragmentation models specify how the energy of the heavy quark is 
transferred to a heavy hadron within the hadronisation process. 

All recent measurements studied the consistency of various models for b quark 
fragmentation to B hadrons [256-258]. The best description of the observed x spectra 
of B hadrons is obtained with the Bowler and with the Lund symmetric models 
implemented within the JETSET Monte Carlo event generator [256,258,259] (see the 
references for details on the fragmentation models). In particular the commonly used 
Peterson heavy quark fragmentation model is found to be disfavoured by the data. 

This result is confirmed in a theoretical study [260] which unfolds the perturbative 
and the non-perturbative part of the description of the x spectrum of B hadrons using the 
Mellin transform technique. Essentially, moments of the distribution are calculated and 
the unfolding is performed in moment space where it simplifies to a simple multiplication. 
The perturbative component is taken from JETSET or from a NLL perturbative QCD 
calculation. In both cases the Lund and the Bowler fragmentation functions are found 
to agree reasonably well with the extracted non-perturbative component while other 
fragmentation functions are disfavoured. 

The radiation of soft gluons from heavy quarks with mass m and energy E is 
predicted to be suppressed within a cone around the heavy quark flight direction of 
opening angle a = m/E ii E » M » Aqcd- Assuming direct correspondence of 
parton and hadron distributions (LPHD as discussed above) allows to make predictions 
in MLLA QCD e.g. for the difference in charged particle multiplicities of light and heavy 
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quark events, 5^ [45,235,261]. The MLLA QCD prediction is independent of cms energy 
and its latest value is 5bi = 4.4 dz 0.4 [261]. In contrast, the so-called naive prediction 
for 5bi considers only the reduction in phase space due to the heavy quark mass. Thus 
the observable 5^ becomes energy dependent and decreases with ^/s. 

Figure EDI (right) shows a summary of measurements of 5u at various cms 
energies [262]. The bands indicate the predictions from perturbative QCD and the 
naive model. The predictions become significantly different at high energies. The recent 
measurements using LEP 2 data by DELPHI [263] and OPAL allow to discriminate 
between the predictions and the naive model is clearly disfavoured by the data. The 
dash-dotted line presents an average of all measurements assuming energy independence: 
(5bi = 3.05 ± 0.19 which is significantly lower than the MLLA QCD prediction. This 
apparent discrepancy is of the order of important corrections beyond the MLLA [261]. 

5.5. Gluon splitting into heavy quarks 

The gluon may split into a quark-antiquark pair; the process g — > qq is one of the 
fundamental processes of QCD. When the virtual mass of the gluon is sufficiently large 
the splitting process can produce a pair of heavy (charm or bottom) quarks. In e+e~ 
annihilation this process is rare, because it occurs only at 0{a^). Neglecting quark 
masses one expects as a rough estimate at ^/s = mzo a rate per hadronic Z° decay 
dqq ~ 1 % where q = b,c. 

The study of gluon splitting to heavy quarks is interesting for several reasons. The 
quark masses provide a natural infrared cutoff in a QCD calculation thus avoiding a 
usually present source of divergence and allowing for reliable predictions. The rates of 
gluon splitting into bb, g^-^;, and into cc, ^^cc, are an important source of uncertainty in 
the measurement of the partial widths for Z° bb or Z° cc decays. In particular 
the Z° —>■ bb vertex is sensitive to possible effects of new physics coupling to the b- 
quark [264] . At high energy hadron colliders a large fraction of the observed B hadrons 
are thought to stem from gluon radiation via g — >■ bb. Such processes are a major 
background to many searches for new particles decaying to heavy quarks like e.g. the 
Standard Model Higgs boson. In order to have precise predictions for the production 
rates of b-flavoured jets from QCD processes at hadron colliders the value of g-^-^^ must 
be known well. The most complete calculation [265] includes resummation of soft 
gluon contributions and predicts gcc = 2.0% and g^y; = 0.18% for as{mzo) = 0.118, 
nic = 1.2 GeV and nib = 5.0 GeV. 

For the measurement of gcc the ALEPH collaboration [266] classifies hadronic Z" 
decays tagged by the presence of a D* meson using the observable AMh = Mh — Ml, 
the difference between the scaled invariant masses of the heavy [Mh) and light {Ml) 
hemispheres of an event. It is expected that events containing the process g cc tend 
to have a larger AM^ due to the necessary presence of an energetic gluon jet. Figure 1^ 
shows the AMh spectrum of events where the D* was found in the heavy hemisphere 
after subtraction of the measured contribution from primary charm production. There 



CONTENTS 



56 



is clear evidence for the process g — > cc consistent with expectations. The other analyses 
from L3 [267] and OPAL [268] use the JADE jet clustering algorithm to find 3-jet events. 
The gluon jet is identified as the lowest energy jet or alternatively by OPAL also as the 
jet which splits into two at the largest value of Hcnt using the JADE algorithm. Events 
with a charm tag in the gluon jet are counted as g — > cc. The results are summarised 
in table H 



The measurements of of ALEPH [269], OPAL [270] and SLD [271] select 4-jet 
events in hadronic Z° decays using the Durham algorithm; between about 10% to 15% 
of the events are selected. The jets are searched for the presence of displaced secondary 
vertices indicating B-hadron decays. Events with two such b-tags in two jets lying closely 
together in angle are considered as candidates originating from g bb. The DELPHI 
analysis [272] measures the rate of bbbb events in hadronic Z° decays and the value of g^^ 
is extracted using a theoretical calculation of BR{7P qqg, g hh) / BR{7P bbbb). 
All results are shown in table IHl 

The averages are calculated based on uncorrected and correlated systematic 
errors [273] extracted from the references given in table IHl The off-diagonal elements 
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Figure 21. Distribution of AMh of events with a D* in the heavy hemisphere after 
subtraction of the measured contribution from primary quarks [266]. The histogram 
shows the Monte Carlo prediction for the process g cc. 
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Table 9. Results from various experiments for ^cc and gj^-g with statistical errors 
and the uncorrelated and correlated systematic uncertainties as extracted from the 
references. The weights found in calculating the averages are also shown. 





9cc 
•10-2 


±stat. 


±syst. 
(uncorr.) 


±syst. 
(corr.) 


weight 


ALEPH [266] 


3.23 


±0.48 


±0.36 


±0.39 


0.12 


L3 [267] 


2.45 


±0.29 


±0.32 


±0.43 


0.20 


OPAL [268] 


3.20 


±0.21 


±0.18 


±0.35 


0.68 


average 


3.05 


±0.16 


±0.14 


±0.36 






5bb 


±stat. 


±syst. 


±syst. 


weight 




•10-3 




(uncorr.) 


(corr.) 




ALEPH [269] 


2.77 


±0.42 


±0.28 


±0.49 


0.38 


DELPHI [272] 


3.30 


±1.00 


±0.41 


±0.67 


0.06 


OPAL [270] 


3.07 


±0.53 


±0.41 


±0.88 


0.09 


SLD [271] 


2.44 


±0.59 


±0.13 


±0.32 


0.47 


average 


2.67 


±0.33 


±0.13 


±0.38 





of the covariance matrix V for the average are calculated from the correlated errors o"i 
and (Tj of experiments i and j as Vy^i^j = min(crj, o"j)^. The resulting weights of the 
individual measurements in the averages are shown in the last column of table IHl The 
results shown here are consistent with those given in [274] but have slightly larger errors. 
The averages gcc = 3.05 ± 0.42% and g^-^; = 0.267 ± 0.052% are larger by two to three 
standard deviations than the theoretical predictions given above. 

6. Inclusive observables 

Inclusive observables are independent of the topology of the hadronic events. Prominent 
examples are -Re+e- = o'{e'^e~ hadrons)/cr(e"''e- — > /i+/i-), the analogous observable 
for Z-decays Rz = r(Z° —>■ hadrons)/r(Z'' — >■ leptons) or the hadronic branching ratio 
in r-decays Rr = r(r hadrons)/r(r leptons). Inclusive observables depend 
on only one energy scale and are effectively QCD corrections to processes predicted 
by the electroweak theory involving quarks. The three observables Re+e-, Rz and Rr 
describe the QCD corrections to inclusive hadronic final states produced in decays of all 
electroweak gauge bosons 7, and W^. A comparison of the QCD corrections for these 
three different electroweak gauge boson decays to hadrons is an important consistency 
test of the theory of strong interactions. 

Since QCD corrections to electroweak decays to hadrons generally scale in LO 
like (1 + as/vr) at high energies careful measurements with uncertainties below (9(1%) 
are needed. The LEP experiments now provide precise measurements of electroweak 
observables such as cross sections from the line shape of the resonance or the 
decay widths of the Z^ boson. Another important set of results from LEP are precise 
measurements of hadronic and leptonic r decays which form the basis for detailed studies 
of the large QCD corrections to the electroweak r decay. Data for precision studies of 
Re+e- come from the low energy experiments [275,276]. 
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6.1. Z° properties 

The LEP experiments have pubhshed their final results on electroweak observables and a 
consistent combination of the data is available [264]. The line shape of the TP resonance 
assuming lepton universality is described by the mass of the Z° boson mzo, its total 
decay width Fz, the total hadronic cross section (Xh, the ratio of hadronic and leptonic 
decay width of the TP boson Ri and the leptonic pole forward-backward asymmetry 
Apg. In addition the partial decay widths to leptons, hadrons and invisible particles 
(neutrinos) F^, F^ and Fjnv are measured. 

The results of the combination are preliminary but since the input values and 
the combination procedure [277] are final the results are not expected to change. The 
important QCD corrections to the electroweak observables derived from the T^ resonance 
are known in NNLO QCD and rather precise determinations of the strong coupling 
asirriio) become possible. 

We concentrate here on a derived set of four electroweak observables with high 
sensitivity to QCD corrections [264]. These are Fh, Rz and the total cross sections for 
lepton and hadron production in e"^e~ annihilation ai and cxh on the peak of the Z° 
resonance. These observables are related to mzo and the partial and total decay widths: 

Fh 127r F,Fh , 127r Vj 

Rz = Tr, = ■ -j;r- and a, = ■ -2 . (46) 

Results for the partial decay widths are known with correlations [264] while correlations 
to mzo are neglected. Using this information we derive the correlations between Fh, Rz, 
(Th and ai as shown in table ITUl 



Table 10. Electroweak observables with high sensitivity to QCD corrections. The 
correlations were derived from correlations for the partial decay width of the Z" 
assuming lepton universality [264] . 









correlations 


Observable 


Result 


Th 


Rz 




rh[MeV] 


1744.4±1.5 


1.0 






Rz 


20.767±0.025 


0.62 


1.0 




crh[nb] 


41.540±0.037 


0.24 


0.18 


1.0 


(T£[nb] 


2.0003±0.0027 


-0.40 


-0.77 


0.48 1.0 



In order to extract values of as{mzo) from measurements of the electroweak 
observables we compare with the complete predictions including electroweak and QCD 
radiative corrections calculated with the electroweak library ZFITTER 6.41 [278]. The 
important QCD corrections are known to 0{al), e.g. for Fh [279,280]. The value of 
as(mzo) is varied between 0.091 and 0.14 and the resulting predictions are interpolated 
until the experimental value is matchedf. The mass of the Higgs boson is taken to 
be mil = 114^45 GeV [264]. Theoretical uncertainties are estimated by a consistent 
variation of the renormalisation scale [281] in the range 0.5 < < 2. In addition mn is 

fThe other input values were nit = 178 GeV and mzo — 91.1875 GeV. 
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varied within its uncertainties. The four results for as(^zo) are shown in table ITT) with 
experimental, theoretical and mn uncertainties and their experimental correlations. 

The theoretical uncertainties for Fh are consistent with an estimate based on 
the difference between the standard approach and an improved theoretical treatment 
which takes so-called vr^-terms into account [282]. An earlier estimate of the 
theoretical uncertainties when using Ri is also consistent with our observed theoretical 
uncertainty [283]. 

A combined value of asirriio) is derived from the individual values as a weighted 
average using the experimental covariance matrix with the theoretical and mn errors 
added to the diagonal. The experimental error of the combined result is obtained 
from the experimental covariance matrix while the theoretical and mn uncertainties 
are found by varying simultanously the renormalisation scale or mn and repeating 
the combination. The combined result is shown in the last row of table ^2 The 
individual values of asirrtio) are consistent within their experimental uncertainties with 
the average and the combination yields a value of x^/d.o.f. = 1.1. The weights of 
the combination are 0.29 for Fh, 0.25 for i?z, 0.17 for dh and 0.29 for ag. Changing 
the combination procedure to use only the experimental covariance matrix or including 
fully correlated theoretical and mn uncertainties also in the off-diagonal elements of the 
covariance matrix yields deviations of the average smaller than 1%. The average value 
is in agreement with the final result for a^{mio) from jet and event shape observables 
discussed in section H] and also with current world averages. 

Table 11. Results for Q;s(?Ti2o) from electroweak observables with high sensitivity 
to QCD corrections. In each row the first error is experimental, the second is from 
variation of and the third is from variation of mpj- The measurements are correlated 
as shown in table EH The average value was obtained as described in the text. 



Observable 


as("T-zo)i exp. ± scale ± mn 


Th 


0.1221 ±0.0037 ±0.0020 ±0.0015 




0.1231 ±0.0037 ±0.0013 ±0.0005 




0.1075 ±0.0069 ±0.0006 ±0.0001 


<yt 


0.1187 ± 0.0030 ± 0.0011 ± 0.0004 


average 


0.1189 ±0.0027±0.0013±0.0007 



This result may be viewed as one of the most reliable measurements of 
as(mzo), because the theoretical uncertainties are small due to the underlying NNLO 
QCD calculations and the inclusive nature of the observables. Uncertainties from 
measurements of the luminosity used to obtain cross sections are included in the 
experimental uncertainties. Uncertainties from non-perturbative contributions are 
expected to be negligible compared to the experimental errors as e.g. the leading power 
correction to the closely related observable Re+e- was shown to scale like where 
Q = mzo is the hard scale of the process [78]. The uncertainty due to the unknown mn 
may be as large as Aas(mzo) ~ 0.003 if ranges of mn up to 1 TeV are considered [207]. 
However, such large Higgs boson masses would have consequences for the formulation of 
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the electroweak theory and the whole procedure of extracting a precise value of a^{mio) 
would have to be reevaluated. 

6. 2. Decay of the r lepton 

Hadronic decays of the r lepton offer unique possibilities to study the physics of strong 
interactions at low energy scales given by the mass of the r lepton rrir = 1.777 GeV. The 
interplay between perturbative and non-perturbative QCD effects is expected to play 
an important role. The ratio of decay widths of the r lepton to hadronic and purely 
leptonic final states Rr = r(r — > hadronsz/^)/r(r — ^ iu^) is an observable like -Re+e- 
and Rz- Leading electroweak effects cancel in the ratio and only QCD and rather small 
electroweak corrections remain. 

A particularity of r decays is the presence of a W boson in the weak decay instead of 
a Z° or 7 in the case of Rz or -Re+e- • This leads to the possibility of hadronic final states 
with even or odd number of pions connected to vector (V) or axial-vector (A) currents 
which has interesting consequences for the size of non-perturbative contributions as 
discussed below. The contributions from perturbative and non-perturbative processes 
may be separated by studying the spectrum of invariant masses of the hadronic final 
states [284]. 

The LEP experiments ALEPH and OPAL have analysed hadronic final states from 
r decays using large samples of Z° decays produced on resonance [285,286]. In the 
ALEPH (OPAL) analysis 124k (150k) r pair candidate events are used. Figure |221 
(upper row) shows the spectral functions v{s) = dRr^v/ds and a{s) = dRr^a/ds with 
RT,v/a{s) = T{t pions Ur) /r(T eVe^r) of hadrouic r decays into final states with an 
even or odd number of pions. The variable s denotes the invariant mass of the hadronic 
final state. In the figures the small corrections based on Monte Carlo simulation for 
unobserved final states are also indicated. In the lower row the sum and the difference 
of the vector and axial vector spectral functions v{s) and a{s) are shown. 

The dashed and solid lines in figure |221 show the predictions by the naive parton 
model and perturbative QCD with massless partons and Q;s(mzo) = 0.122. The pQCD 
predictions are seen to agree well with the inclusive v{s) + a{s) spectral function while 
for v{s) and a{s) the prediction lies below or above the data. This observation already 
indicates that there are significant non-perturbative contributions to v{s) and a{s) 
which, however, largely cancel in the sum v{s) + a{s). 

For the detailed analysis moments of the spectral functions are determined [284]: 




for kl = 00, 10, 11, 12, 13. These moments are constructed such that they are 
theoretically and experimentally well under control. The moment R^^^^^{sq = m^) = R^ 
corresponds to the usual hadronic branching ratio in r decays. 
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Figure 22. (upper row) Spectral functions for hadronic r decays to final states 
with even (v(s)) or odd (a(s)) number of pions. The shaded areas show contributions 
from hadronic final states as indicated and the Monte Carlo (MC) based correction for 
unobserved final states, (lower row) Sum or difference of the spectral functions. The 
lines show predictions as indicated on the figures. All figures from [286]. 



The theoretical prediction is formulated as a sum of perturbative and non- 
perturbative QCD contributions: 

R^\v/a{So) = -SEw\Vud\^ (l + Sew + "^pert + ^non-pcrt,»;/a) (48) 

where Sew = 1-0194, Kd = 0.9745 ± 0.0004 [18] and = 0.001. The term 

5pgrt contains the perturbative QCD part known to C(ag) [27,287,288]. The original 
calculation of [288] has been improved in [289] by a more accurate treatment of the 
running strong coupling in the calculation using the full /3-function instead of an 
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approximation. A third calculation includes a resummation of the effects of renormalon 
chains [290]. The non-perturbative contribution can be analysed in the framework of 
the Operator Product Expansion (OPE) as a power series in terms of [288] 

"non-pert, ti/a / ^ ^v/a ' 

D=2,4,6,... 

The first term [D = 2) in equation (j^9|) is a correction for finite parton masses and 
can be calculated perturbatively. The D = 4,Q,8 terms are studied as non-perturbative 
corrections in the analyses. 

As already indicated by figure |221 the non-perturbative corrections turn out to be 
consistent with zero when moments of the inclusive spectral function v{s) + a{s) are 
analysed. The correction terms are found to be of roughly the same size but opposite 
sign when moments of v{s) and a{s) are studied individually [285,286]. 

The analysis of the perturbative term 5^'^^ allows a precise determination of 
the strong coupling constant at energy scales given by the mass of the r lepton 
rrir = 1777.0 ± 0.3 GeV [18] essentially free of non-perturbative uncertainties. 

Following [285,286] the determination of as from hadronic r decays using Rj^ — R 
is updatedf for recent values of the r lifetime t^- = 290.6 ± 1.1 ps and the branching 
ratio = B{t ^iv^Vt) = 0.1736 ± 0.0006 [18]. The contribution of strange hadronic 
T decays is subtracted using Rr^s = 0.1677 ± 0.005 [291]. We find for the non-strange 
branching ratio Rr = 3.468 ± 0.011. A similar study [292] using preliminary ALEPH 
data [293] finds Rr = 3.471 ± 0.011 in agreement with our result. The improved 
calculation of [289] is chosen to extract a value for as- The 0(ag) term is set to 
zero, because for consistency with other determinations of as discussed in this report 
we prefer to extract a value of as without the additional assumptions needed to estimate 
the O(a^) term [294] f. We find 

as(m^) = 0.347 ± 0.005(exp.) ± 0.003(non - pert.) ± 0.018(theo.) (50) 

The non-perturbative correction was set to 5non-pert,»;+a = —0.0024 ± 0.0025 [286] 
(consistent with 5non-pert,i)+a = —0.0048 ± 0.0017 [293]) and varied within its errors 
in order to obtain the uncertainty due to non-perturbative effects. The theoretical error 
was determined by a variation of the renormalisation scale in the range 0.5 < < 2.0. 
We take the larger variation from changing to define the theoretical uncertainty. We 
did not consider a variation of the renormalisation scheme as done in [286] in order to 
be consistent with other determinations of as discussed in this report. After evolving 
to the Z° mass scale using the 3-loop /^-function and 2-loop matching at the crossing of 
heavy quark thresholds [20] we find 

ag(mzo) = 0.1221 ± 0.0006(exp.) ± 0.0004(non - pert.) (51) 
± 0.0019(theo.) ± 0.0003(evol.) . 

fCode and help by S. Menke is greatly appreciated. 

|In [294] the rij and contributions to the 0{a'^) term are given. Further parts of the 0{a'^) 
term are estimated using the renormalisation scheme optimisation techniques PMS and ECH (see 
section 113^ . 
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The last error (evol.) is due to uncertainties in the evolution procedure. The evolution 
of as(?7T'r) to aQ,{mio) takes the crossing of the charm and bottom quark threshold into 
account by setting the charm quark mass to rUc = 1.778 GeV. The uncertainty due to 
the evolution procedure is evaluated following [295] by setting the matching thresholds 
to mc = 4 GeV or rub = 20 GeV and by changing the absolute quark masses within 
their uncertainties [18]. 

The analogous determination of a;s(^r) and as{mzo) using the calculation of [290] 
results in as(mr) = 0.303 ± 0.003(exp.) ± 0.002(non - pert.) ± 0.011(theo.) and 
as(mzo) = 0.1166 ± 0.0004(exp.) ± 0.0003(non - pert.) ± 0.0015(theo.) ± 0.0003(evol.). 
The difference between the result in equation (jKT|) and this result is Aas(mzo) = 0.0055 
which corresponds to less than three standard deviations of the total errors. However, 
since the calculation of [290] is rather different compared to e.g. [289] such discrepancies 
might be expected. In [214] the strong coupling constant is extracted using the ALEPH 
data [285] and a renormalisation group invariant formulation of the perturbative QCD 
prediction. Their result is as (^zo) = 0.1184±0.0007(exp.)±0.0019(theo.)±0.0006(evoL) 
which is in fair agreement with ()51|) . A similar result was found in [296]. The 
difference between using the calculation of [288] and the improved calculation of [289] 
is Aa;s(mzo) ~ 0.003 [286]. Using the prediction of [289] including the estimate of the 
0{al) term as in [286] we find Aas(^zo) = -0.0006. 

This measurement of as{mzo) shown in equation 1)511] has comparable uncertainties 
to and is consistent with the measurement using the Z° line shape observables discussed 
in section 16.11 Even though the difference between different methods of extracting 
as{fnzo) from the hadronic r data appears large it corresponds to about three times 
the total uncertainties. We conclude that the determination of as{mzo) from hadronic 
r decays is reliable within the reasonably well estimated total errors. 

6.3. e~^e~ annihilation at low energies 

In the e"'"e~ annihilation to hadrons at low energies the intermediate gauge boson is 
predominantly a photon thus giving access to the process 7 — > qq. The quantity of 
interest is -Re+e- = (j{e'^e~ hadrons)/cr(e"'"e^ fi^fi^), because many experimental 
and theoretical uncertainties introduced by selection efficiency corrections or radiative 
corrections are reduced in the ratio. Many measurements exist and have been 
summarised in other reviews, see e.g. [18,197,297]. Theoretical predictions for this 
observable exist in C(ag) [27,36]. 

In [197] a consistent analysis of -Re+e- using TRISTAN and PETRA data [298] is 
used to extract the strong coupling constant in 0{al) with the result as (42.4 GeV) = 
0.175±0.028 or as(mzo) = 0.126±0.022. The measurements of i?e+e- take account of the 
Z° mass measured at LEP and also include improved electroweak radiative corrections. 

The CLEO collaboration has measured -Re+e- at y/s = 10.52 GeV, i.e. just below 
the Y(4S) resonance. A QCD prediction including heavy quark masses and corrected 
for QED radiation is used to extract as(10.52 GeV) = 0.20 ± 0.01(stat.) ± 0.06(syst.). 
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For the evolution to the Z° mass scale we use a b quark threshold of 10.53 GeV such 
that the value as (10.52 GeV) = 0.20 is interpreted as the strong coupling in a regime 
with four active quark flavours. The result is 

as(mzo) = 0.130 ± 0.004(stat.) ± 0.025(syst.) (52) 

in good agreement with previous results. The dependence of this result on the actual 
value of the b quark threshold and the b quark mass = 4.25 ± 0.15 GeV [18] in the 
MS scheme is negligible. 

The behaviour of R^+e- {s) at low values of s in the region of hadronic resonances 
may be described by perturbative QCD in the same way as the spectral functions from 
hadronic r decays [275]. The perturbative prediction for R^+e- is written in terms of a 
the so-called Adler D-function: 

1 r D(s) 

i?e+e-(so) = — / -^ds (53) 
2m Js=so s 

where the D-function D{s) is known in 0(as) [27,36]. Non-perturbative effects are 
treated by applying the OPE to moments of Re+^-^s): 




Similar to the case of hadronic r decays the complete prediction for D{s) is a sum of 
perturbative and non-perturbative termsf. The non-perturbative part of D{s) in the 
OPE is a power series in terms of 1/s" starting at n = 3. 

In [276] the same data as in [275] are used to determine the moments 
(4 GeV^), kl = 20, 30, 31, 32, 33. A simuhaneous fit of the QCD prediction for these 
moments with as (2 GeV), the strange quark mass and the first two non-perturbative 
terms proportional to 1/s^ and gives a good description of the data. The result 
for the strong coupling constant is as(2 GeV) = 0.286 ± 0.031(exp.) ± 0.015(theo.) 
corresponding to 

as(mzo) = 0.117 ± 0.005(exp.) ± 0.002(theo.) . (55) 

Compared to the result from hadronic r decays shown in equation ()51|) the theoretical 
uncertainties are similar while the experimental errors are larger. Within the 
uncertainties the measurements from hadronic r and decays and from hadron 
production in e"'"e~ annihilation are consistent. New measurements of hadron production 
in e"'"e~ annihilation at low ^/s promise to allow reduced experimental uncertainties in 
the future [299] such that this measurement should have competitive uncertainties. 

7. Summary of as determinations 

In this section we collect and systematically compare all significant measurements of the 
strong couphng constant as(mzo) discussed in this report. We concentrate on reliable 

f Corrections for finite quark masses are viewed as part of the perturbative prediction. 
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measurements with comparatively small uncertainties in order to probe the consistency 
of the theory in the best possible way. 

Table [T21 presents the results selected from various sections of this report. The 
first two rows of the table contain the results based on C(ag)+NLLA QCD with 
hadronisation correction performed using Monte Carlo simulation. As discussed in 
section 14.4.11 the analyses from JADE and LEP use the same set of 3-jet event shape 
observables sensitive to as in 0{as)- The ALEPH, DELPHI and OPAL analysis (see 
section I4.4.2|l based on the 4-jet fraction R4 are sensitive to as in 0{a^). ALEPH 
and OPAL use NLO+NLLA calculations and the Durham algorithm while DELPHI 
uses NLO calculations with experimentally optimised renormalisation scale and the 
Cambridge algorithm. 

A direct comparison of the results from JADE and LEP shows very good agreement. 
As the set of observables is identical and the analysis procedures are rather similar the 
good agreement shows that a consistent analysis of 3-jet event shape observables is 
possible from the lowest PETRA/JADE cms energies of 14 GeV to the highest LEP 2 
cms energies of about 200 GeV. The results from analysis of R4 are consistent with 
the other NLO values within the uncertainties. This confirms that the calculation of 
radiative corrections within perturbative QCD is consistent within this set of rather 
different observables. 

Table 12. Summary of determinations of as("izo) based on observables, theoretical 
calculations and experiments as indicated. The experimental errors (exp.) also include 
statistical errors. Numbers covering two columns represent combined errors. The 
observables R4(D) and R4(C) refer to R4 using the Durham or Cambridge jet finding 
algorithm. 







Experiments 








Observable 


Theory 


(Q [GeV]) 


as(mzo) 


±exp. 


±had. ±theo. 


3-jet shp. 


NLO+NLLA 


JADE (14-44) 


0.1203 


0.0018 


0.0053 0.0050 


3-jet shp. 


NLO+NLLA 


LEP (91-207) 


0.1201 


0.0009 


0.0019 0.0049 


R4(D) 


NLO+NLLA 


ALEPH (91) 


0.1170 


0.0009 


0.0010 0.0017 


R4(D) 


NLO+NLLA 


OPAL (91-207) 


0.1182 


0.0015 


0.0011 0.0018 


R4(C) 


NLO 


DELPHI (91) 


0.1175 


0.0011 


0.0027 0.0007 


sc. viol. 


NLO 


LEP (14-207) 


0.1192 


0.0056 


0.0070 


C^L/c^tot 


NLO 


LEP (91) 


0.1169 


0.0035 


0.0018 0.0072 


Th, -Rz, 0"h, 0"£ 


NNLO 


LEP/SLD (91) 


0.1189 


0.0027 


0.0015 




NNLO 


e+e" (1.777) 


0.1221 


0.0006 


0.0004 0.0019 


-Re+e- 


NNLO 


e+e- (2) 


0.117 


0.005 


0.002 



Having observed good consistency between the 3-jet results we construct an average 
value of the results from JADE and LEP assuming that statistical errors are uncorrelated 
and experimental and hadronisation errors are partially correlated. Theoretical errors 
are treated as in the combinations of as in section 14.4.11 The result for as(mzo) from 
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3-jet event shapes using NLO+NLLA QCD predictions is 

ag(m2o) = 0.1202 ± 0.0004(stat.) ± 0.0009(exp.) (56) 
± 0.0025(had.) ± 0.0049(theo.) . 

The weights are 0.32 for the JADE result and 0.68 for the LEP resuh. The determination 
of as(wT.zo) in 0(ag)+NLLA from 3-jet event shape observables has a fairly large 
theoretical uncertainty compared with some other measurements. This uncertainty 
might be reduced substantially when complete NNLO calculations for 3-jet production 
in e"'"e~ annihilation will become available [25,26]. 

The measurements of a^{mio) from R4 discussed in section 14.4.21 have small 
uncertainties. However, since no other comparable measurements based on different 4-jet 
observables exist consistency checks were not possible. For the purpose of presentation 
we combine the ALEPH and OPAL measurements based on R4(D) using the Durham 
algorithm and NLO+NLLA theory and the DELPHI measurement based on R4(C) 
using the Cambridge algorithm and NLO theory. Experimental errors are assumed 
as partially correlated while hadronisation and theoretical uncertainties are treated by 
simultaneous variation of the input values within their individual errors. The result from 
analyses of 4-jet fractions R4 becomes 0:3(^20) = 0.1175±0.0002(stat.) ±0.0010(exp.) ± 
0.0014(had.) ± 0.0015(theo.). 

The result from the NLO analyses of scaling violation (sc. viol.) in e+e~ annihilation 
to hadrons is discussed in section 15.21 The method is inclusive and thus has negligible 
hadronisation uncertainties. There is good agreement with the results from other 
observables. However, since the QCD calculations are limited to NLO a direct 
comparison and combination with other inclusive measurements of 0:3 (mzo) based on 
NNLO QCD discussed below is not attempted. 

The seventh row of table ^1 gives the result of the NLO analysis of the ratio of 
longitudinal to total cross section crL/fXtot from section Combining all available 
measurements of this observable and using the hadronisation corrections derived 
from [247] yields a consistent measurement of Q!3(mzo) with a significantly reduced 
experimental error compared with [248]. This measurement is combined with the 
determination of aQ,{mio) from scaling violations assuming experimental errors as 
partially correlated while hadronisation and theory errors are treated by simultaneous 
variation of the input values within their hadronisation or theoretical uncertainties. We 
find for a^[m'zp) from analysis of fragmentation in e"'"e~ annihilation to hadrons: 

ag(mzo) = 0.1179 ± 0.0040(exp.) ± 0.0010(had.) ± 0.0071(theo.) . (57) 

The lower three rows of table ^1 present results from the inclusive observables 
discussed in section IHl derived using 0{a\) (NNLO) QCD calculations. The inclusive 
nature of the observables allows the calculation in 0{a^ containing two-loop radiative 
corrections and reduces uncertainties from non-perturbative effects. The results for 
03(77120) from the Z° line shape observables Fh, -Rz, CTh and at and from analysis of Rr 
are in agreement within less than two standard deviations of their total uncertainties. 
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The good agreement of the precise values of as('^zo) from the TP hne shape observables 
and Rr at the level of Aas("^zo) = 0.0032 ± 0.0032 may be viewed as a strong test of 
asymptotic freedom, since two measurements obtained at energy scales differing by a 
factor of (9(100) are compared. In the calculation of the error of Aas(^zo) only the 
smaller of the theoretical uncertainties was considered. 

We derive an average value of the NNLO results assuming statistical errors as 
uncorrelated and experimental errors to be partially correlated. Hadronisation and 
theoretical uncertainties are treated by simultaneous variation of the input values within 
their uncertainties. The result for Q;s(mzo) from inclusive observables in NNLO QCD is 

ag(mzo) = 0.1211 ± 0.0010(exp.) ± 0.0018(theo.) (58) 

where statistical errors have been included in the experimental error and hadronisation 
uncertainties are included in the theory error. The weights are 0.27 for the TP line shape 
observables, 0.70 for Rr and 0.03 for -Re+e- with x^/d.o.f. = 1.05/2. The consistency 
between the three results is impressive as indicated by the x^/d.o.f. value. This result is 
considered to be the most reliable determination of the strong coupling constant ag(mzo) 
discussed in this report, because it is based on NNLO QCD calculations and has only 
small uncertainties from modelling of non-perturbative effects. Further improvements 
of this determination of asijnio) may be expected from new measurements of -Re+e- in 
the low energy region leading to a reduction of the experimental uncertainties [299]. 

Figure ESI summarises the high precision measurements of as discussed in this 
section. The data are compared with our average of NNLO QCD based determinations 
of as(^zo) given in equation (|58|l shown as lines and grey band. The value of as(^zo) is 
evolved in NNLO to other cms energies as explained in section 12.11 The consistency 
between the predicted behaviour of as(yi) using only the NNLO observables with 
the other measurements is good; this shows that QCD successfully describes many 
important features of hadron production in e"'"e~ annihilationf at scales covering two 
orders of magnitude from rrir to the highest LEP 2 energies of about 200 GeV. The 
most significant difference is Aas('^zo) = 0.0035 ± 0.0021 between the combined results 
from R4 and from inclusive observables in NNLO. Both errors are combined taking 
experimental and hadronisation uncertainties as uncorrelated and theoretical errors as 
fully correlated. Since R4 is an exclusive observable based on 4-jet final states in contrast 
to the fully inclusive observables assuming experimental and hadronisation uncertainties 
to be uncorrelated seems reasonable. 

8. Quark and gluon jets 

Differences between the properties of jets originating from quarks or gluons are a basic 
consequence of QCD. The quarks and gluons can be viewed as carriers of colour charges 
given by their colour factors which are Cp = 4/3 for quarks and Ca = 3 for gluons. 

f More generally speaking, the analyses cover hadronic decays of electroweak gauge bosons. 
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Figure 23. Values of as determined at cms energy from various processes as 
indicated. The NNLO incl. data are from Rr{mr), i?e+e-(2 GeV) and the Z° line 
shape observables. The error bars give experimental and total uncertainties. Invisible 
error bars are smaller than the symbols. The measurements at ~ mzo have been 
separated horizontally for clarity. The lines and shaded band show the running of the 
strong coupling based on the average value shown in equation l|58|l . 



Predictions for e.g. the ratio of parton multiplicities 

(59) 



_ Ca 



in), Cp 

in jets originating from hard gluons or quarks [300, 301] can be related to observed 
quantities like charged hadron multiplicities measured in experimentally identified quark 
or gluon jets by assuming local parton hadron duality [64]. One therefore expects the 
particle multiplicities in gluon jets to be larger than in comparable quark jets. 

With the relations {x)^{n)^ = l,p = q,g, where x = Pparton/Pjet is the scaled 
momentum of a parton produced in a quark or gluon jet and (n)^ is the parton 
multiphcity in the jet, one gets {x)^/{x)g = Ca/Cf- Thus gluon jets can be expected 
to have a softer scaled hadron momentum spectrum (fragmentation function) compared 
to quark jets of the same energy [301]. The tails at large a; ^ 1 of the fragmentation 
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functions in gluon or quark jets can be predicted to have a ratio ~ (1— x)/ln(l— [301]. 
This imphes that the gluon jet fragmentation function is expected to be softer by a factor 
of approximately 1 — x. 

In experimental studies of differences between jets stemming from hard quarks or 
gluons the main difficulty lies in achieving a clear and unambiguous identification of the 
original parton. Due to the confinement property of QCD free partons are not observed 
but instead they form more or less collimated showers of softer partons and in turn 
of hadrons. This may cause a loss of the close correspondence between the hadronic 
final state and the underlying hard parton state. In addition the correlated colour fiow 
between the quarks and gluons leads to effects of colour coherence with the consequence 
that the properties of gluon jets depend on the topology of the state in which they were 
produced, see e.g. [302]. 

Most theoretical predictions of differences between e.g. the average multiplicity 
of charged hadrons (nch.) in quark or gluon jets don't take experimental definitions 
of jets into account [303]. The jets are simply defined as qq or gg systems produced 
from a colour singlet point source. Hadronic final states in e"'"e~ annihilation are a 
source of qq systems but gg systems from a colour singlet point source rarely occur 
in nature, the decay Y(3S)— > 7x(2P);,^j=o,2, x(2-P)fe,j=o,2 ^gg [304] or production via 
2-photon interactions [305] being examples. This prevents a clean interpretation of the 
experimental results based on jet finding algorithms to define jets and thus experimental 
techniques matching the theory more closely had to be found. 

The next sections describe results obtained using exclusively defined jets employing 
jet finding algorithms and results based on inclusive jet definitions without involvement 
of jet finding algorithms. 

8.1. Exclusive jets in 3- jet events 

The first analysis showing significant differences between quark and gluon jets used LEP 
data from the TP peak and a novel technique to obtain clean samples of quark and gluon 
jets of the same energies [306]. The method [307] employs lepton tagging of quark jets 
in one-fold symmetric 3-jet events, so-called Y-events. The selection demands the angles 
between the highest energy jet and the other two jets to be within 150° ± 10° yielding jet 
energies of about 24 GeV for the lower energy jets. The presence of a high momentum 
lepton (electron or muon) is required in one of the two lower energy jets. Since most 
such leptons stem from weak heavy quark decays the corresponding jet is identified as a 
quark jet. The production of heavy quark pairs from gluons of the perturbative cascade 
contributes (3.1 ± 0.4)% to c and (0.27 ± 0.05)% to b production and is thus a small 
effect (see section I^TSj) . The most energetic jet is also a quark jet with high probability 
due to the bremsstrahlung nature of hard gluon radiation and thus the remaining jet 
is likely to be a gluon jet. The selection of a second event sample requiring the same 
selection criteria except the lepton tagging allows to correct measurements of quark jet 
properties for a possible lepton tagging bias. The gluon jets are assumed to be unbiased. 
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Figure 24. (left) (a) Scaled energy spectrum of particles from quark (points) or gluon 
(line) jet cores. The solid points are corrected for the tagging bias, (b) Ratio of gluon 
to corrected quark jet data of (a). The lines in (b) are predictions from JETSET 7.2. 
(right) (a) and (c) Multiplicity spectrum as a function of the angle in the event plane 
■0 for quark (points) and gluon (line) jets on a linear or logarithmic scale. The solid 
points are corrected for the tagging bias. The angle ^p runs for points on the quark jet 
side and for the line on the gluon jet side, (b) Gluon/quark jet ratio of the data in (a) 
and (c). All figures from [306]. 



The main results of tlie analysis are summarised in figure EI On the left in (a) the 
scaled energy spectrum of particles in the cores of the quark or gluon jets is shown where 
the open (solid) points represent the data after lepton tagging (after lepton tagging and 
bias correction). The line displays the result for gluon jets. The jet cores are defined as 
the angular region 135° < ip < 165° where ip is the angle in the event plane w.r.t. the 
highest energy jet. In left (b) the ratio of the bias corrected quark to the gluon jet data is 
given. The plots show convincing evidence for a softer scaled energy spectrum in gluon 
jets, as predicted. On the right in (a) and (c) the multiplicity spectrum obtained in 
quark (points) and gluon (line) jets is presented as a function of the angle ip- The open 
and solid points are defined as above. The angle ijj runs on the quark jet side of the event 



CONTENTS 



71 



plane for points and on the gluon jet side for the hne. Figure (b) on the right shows the 
gluon/quark jet ratio for 112.5° < ifj < 180°. The figures show evidence for a broader 
particle distribution in gluon jets compared to quark jets. The discrepancy visible in 
(c) in the region between jets is a manifestation of the string effect, see e.g. [308]. 

In a similar analysis using the same technique of lepton tagging of quark jets in 3-jet 
events the OPAL collaboration studied distributions of electric charge in quark and gluon 
jets [309]. The electric charge of a jet is defined by the sum Q of charges of all charged 
particles assigned to the jet. The distributions of Q filled separately for events with a 
positive or negative lepton tag are observed to have mean values significantly different 
from zero for the tagged and the highest energy jet. The Q distribution for gluon jets has 
a mean value consistent with zero. The signs of the quark jet Q distribution mean values 
correspond to the expected charge of the primary quark or antiquark. This observation 
constitutes experimental evidence based on e"'"e~ data that gluons don't have electric 
charge. 

These analyses have been refined in several ways. The use of vertex instead of 
lepton tagging of heavy quarks employs that the decay vertex of heavy hadrons can 
be reconstructed using micro-vertex detectors. The technique uses a larger fraction of 
heavy hadron decays and allows better efficiencies and purities. The use of alternative jet 
finding algorithms like the cone [186] or the JADE algorithm quantified the influence 
of the jet reconstruction method on the results [310,311], in particular the ratio of 
average charged particle multiplicities in quark and gluon jets at 24 GeV jet energy was 
found to be biased by about 10%. A different correction method (qqg/7) uses hadronic 
events consisting of two jets and one isolated and energetic photon as a pure sample 
of quark jets [311]. This method allows analysis of quark and gluon jets at varying 
energies finding e.g. that the ratio of charged particle multiplicities in gluon and quark 
jets depends on the jet energy. 

A direct comparison of quark jets originating from light (u, d or s) or b quarks 
with gluon jets was performed by tagging the high energy jets of the one-fold symmetric 
samples [73,312]. Vertex tagging is used to select b jets while the absence of charged 
particle tracks with large positive impact parameter indicates light quarks. Comparing 
the gluon jet tagged sample (vertex tag in a lower energy jet) with the light or b tagged 
sample (impact parameter or vertex tag in the high energy jet) allows to extract the 
properties of pure gluon, light or b jets at jet energies of about 24 GeV. Some results 
for the ratios of charged particle multiplicities are shown in table El 

The multiplicity ratios -Rg/q and -Rg/uds are significantly larger than one in all cases 
proving that the charged particle multiplicity is greater in gluon than in quark jets. 
However, the difference between the ratios -Rg/q and -Rg/uds based on the different 
jet finding algorithms (Durham, JADE or Cone) implies that the results cannot be 
compared directly to QCD predictions [73,303]. The values of -Rg/b show that at jet 
energies of ~ 24 GeV the multiplicities of gluon and b jets are comparable. 

Figure ESI (left) [73] shows the differential energy profile 0£;(r/-R) = AE/A{r/R) 
where A-E" is the energy contained in an annulus of width A{r/R) at half angle r aligned 
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Table 13. Ratios of charged particle multiplicities of gluon and inclusive quarks 
(i?g/q), gluon and light quarks (i?g/uds) or gluon and b quark jets (i?g/b) at ~ 24 GeV 
jet energy based on the Durham, JADE or the Cone algorithm. The methods "Y" and 
"qqg/7" refer to the analysis methods based on one-fold symmetric Y-events or qqg 
and qq7 events to obtain results for pure quark and gluon jet samples from the data. 



Jet alg. 


Vent 


method 


refs. 




^g/uds 


Rg/h 


Durham 


0.01 


Y 


[312] 


1.19 ±0.03 


1.25 ±0.09 


1.06 ±0.05 


Durham 


0.02 


Y 


[73,310] 


1.25 ±0.04 


1.39 ±0.05 


1.09 ±0.03 


Durham 


0.015 


Y 


[311] 


1.25 ±0.04 






Durham 


0.01 


qqg/7 


[311] 


1.24 ±0.06 






JADE 


0.04 


qqg/7 


[311] 


1.35 ±0.08 






Cone 




Y 


[73,310] 


1.10 ±0.03 


1.14 ±0.04 


0.92 ±0.04 



with the axis of the Cone jet with half angle R. A broader energy profile of gluon 
compared to light quark jets is clearly visible. The energy profile of b jets at ~ 24 GeV 
jet energy turns out be comparable to the gluon jet profile. Figure EH (right) [73] 
presents the fragmentation function (FF) of charged particles. The gluon FF is clearly 
softer than the light quark (uds) FF, but the b quark FF is again observed to be similar 
to the gluon jet FF at jet energies of about 24 GeV. 




Figure 25. (left) Comparison of b, uds and inclusive quark jets and gluon jets using 
the cone jet algorithm for the differential energy profile 4>E{r/R, see text, (right) 
Comparison of the same b, uds and inclusive quark jets and gluon jets for the charged 
particle fragmentation function. All error bars show statistical uncertainties. Both 
figures from [73]. 



All studies discussed so far did not take into account the topology dependence of 
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jet properties, see e.g. [45], because with symmetric 3-jet events the effects of topology 
dependence are suppressed. The topology dependence of particle multiplicity in jets in 
3-jet configurations has been studied theoretically in detail, see e.g. [313] and references 
therein. 

In the MLLA the total particle multiplicity of 3-jet events A^qqg(s,p_L) produced 
at cms energies a/s also depends on the transverse momentum of the gluon jet, 
due to colour coherence [45,314]. The quantity is expressed as a sum of particle 
multiplicities in biased qq events A^qq(s,p_L) and in ideal and unbiased gg events iVgg(px) 
taken at appropriate scales [313,315]: 

iVqqg(s) = N^Ti{L,P^,Lu) + \Ngg{p^^Lu) (60) 

with L = \n(s/AQQj)). The contribution Nq^{L,pj_) carries an explicit dependence on 
the gluon jet scale p± which relates it to the unbiased particle multiplicity Nq^ measured 
without jet requirements [315]: 

iVqq(L,p^) = N,^{L') + {L- L')^^^^ , (61) 
with L' = 2 ln(p_|_/AQCD) + 3/2). The transverse momentum of the gluon jet is defined 



by p±,Lu = Y-^qg-Sqg/s with Sqg and Sqg the invariant masses in the qg and qg systems. 
In an alternative expression for A^qqg the variable L is replaced by Lqq = ln(sqq/AQQj-,) 



and p_L^Lu is replaced by p_L,Le = ySqgSqg/sqq [45,314]. 

These predictions are only valid for 3-jet events selected without a cut on a fixed 
value of the jet resolution parameter ?/cut, since such a cut introduces an upper limit on 
the transverse momenta of any subjets. In contrast, when three jets are reconstructed 
in all events and the events are classified by the p± of the third jet the only remaining 
bias is given by the limitation of subjet transverse momenta. This bias is taken into 
account by the quantity iVqq discussed above. Also, the jet reconstruction should be 
done using a jet resolution definition based on transverse momentum, e.g. the Durham, 
Cambridge or LUCLUS/PYCLUS algorithms. 

For 3- fold symmetric 3-jet events with all three angles between the jets ~ 120°, so 
called Mercedesf events, the alternative expression discussed above for the transverse 
momentum of the gluon jet simplifies [45,314]: 

This has been used as a motivation in experimental studies to use gluon and quark jet 
scales Qjet = -Ejet sin(6/2) where 6 is the angle to the closest jet and thus takes the 
event topology into account [316-319]. The study [316] showed that using Qjet for quark 
jets and Qje^ = ^QqgQqg removes biases in the jet multiplicities which appear when jets 
with similar energies Ej^t but at different angles B are compared. 

The theoretical description of multiplicity in 3-jet events of equations (jUUj) and (pT|l 
has been tested. In [320] the charged particle multiplicity in 2-jet events selected with 
f Trademark of DaimlcrChrysler AG acknowledged 
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the Durham jet algorithm with fixed vales of ?/cut has been measured. The measurements 
have been compared with equation using the relation = -E'visi/Z/cut and reasonable 
agreement within the total uncertainties of the data has been observed. Based on this 
result the charged particle multiplicity in unbiased gg events has been extracted as a 
function of the jet energy scale Q. The results are shown in figures EEl (a) and (b). In 
figure ESI (a) the measured Nqq{L,p±) have been used to extract Ngg shown by the solid 
points at Q = p±. The measurement of Ngg is compared with a pQCD prediction [321] 
which has been fitted to the data (solid line) and predictions by Monte Carlo simulations 
(dashed and dash-dotted lines) of artificial gg events. The predictions are found to 
agree well with the data including other measurements from inclusively defined gluon 
jets shown by solid triangles (see below). In figure 12^1 (b) the Nqq{J,p±) have been 
calculated using equation (j6T| using both definitions of p±; stars correspond to p±^Lu 
and J = L and open points to p±,Le and J = Lqq. The Monte Carlo predictions and 
inclusive gluon jet data are the same as in figure EHl (a). The results based on p±^Lu are 
in better agreement with the inclusive gluon jet data at lower and higher Q than the 
results employing p±^Le as the jet energy scale. 

The quantities Ngg for gluon jets and Nq^ for quark jets in two-parton systems are 
related by diVgg(L')/dL = ^^/^^(l -aoCr/^)diVqq(L)/dL with ao = 6Ca/ illCA-2nf) 
and Cr = 107r^/27 — 3/2 [315]. The measurements of N^^ together with measured iVqq 
thus allow a precise determination of ratio of colour factors Ca/Cf [320]: 

Ca/Cf = 2.23 ± 0.01(stat.) ± 0.14(syst.) (63) 

The data in figure 1211 (a) above Q = 7 GeV have been used in the fit and the fit 
result corresponds to the solid line. In [318] a similar measurement has been performed 
for the first time using slightly different predictions [235] and involving an additional 
non-perturbative correction. The result was Ca/Cf = 2.25 ± 0.06(stat.) ± 0.12(syst.). 
Both results for Ca/Cf agree well with the QCD expectation of Ca/Cf = 2.25 and 
have small systematic uncertainties compared to other measurements of this quantity. 
A compilation of other measurements of Ca/Cf is shown in sectional 

The FFs of quark and gluon jets in mirror symmetric Y-events and completely 
symmetric Mercedes events were studied as function of the scaled jet energy Qj^t [317]. 
The comparison of the gluon FF in Y and Mercedes events found evidence for scaling 
violations in the gluon FF. The scaling violation of the gluon jet FF was seen to 
be stronger than the comparable quark jet FF consistent with the QCD expectation 
motivated by the larger colour charge of gluons compared to quarks. This observation 
was confirmed by a study of subjets in quark and gluon jets. To find subjets the jet 
finding algorithm is run only on the particles associated with the gluon or quark jet [317]. 
The fraction of events without subjet production was found to be smaller in the quark 
jet sample in agreement with the expectation and an NLLA QCD calculation. 

Other subjects of study have been the production rates of identified particles in 
quark and gluon jets. In several analyses the production of tt^ and rj [319,322], of 
K° [319,323-325], K+ and p [324], and A° [323-325] was investigated. In [322] the 
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Figure 26. (a) Results for iVgg shown by solid points as function of the jet energy scale 
Q based on measured Nqq{L' ,p±) [320]. The error bars indicate total uncertainties, 
(b) Results for A/gg based on calculated Nqq{L',pj_) using (stars) or pj_^Le (open 

points). Both figures show by the solid and dashed or dash-dotted lines the same 
predictions from pQCD and Monte Carlo simulations and by the solid triangles other 
measurements based on inclusive gluon jets [320]. 



production of vr'^ and rj was measured in 2- and 3-jet events and rj production was 
found to increase stronger than production in 3-jet events. This observation was 
not confirmed in [319] where production rates of 7r°, rj and K3 in gluon jets were found 
to be larger than in quark jets by the same factor as for the increase of the inclusive 
charged particle multiplicity. In [324] it was observed that the momentum spectra of Kg, 
K"*", p and A*^ were softer in gluon compared to quark jets with a higher total particle 
multiplicity. 
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8.2. Inclusive jets 

In the definition of inclusive gluon and quark jets no jet finding algorithm is directly 
involved. This removes systematic uncertainties connected with the assignment of 
particles to jets, as discussed in the previous section. Studies using inclusive gluon 
jets are based on i) particle decays to hadronic final states via gluons only [304, 326] or 
ii) e"'"e~ qqg events at high energy in the rare kinematic configuration where a hard 
gluon g recoils against the qq system [303] . 

The first study using inclusive gluon jets was [304] using the decays Y(3S)— >■ 
7X6(2P), Y(1S)— >-ggg, and qq events produced in the continuum at y/s = 10.55 GeV. 
The Xfe(2P) states with J = 0, 2 can decay to gg while the Xfe(2P), J = 1 state decays 
to qqg. The three different Xf,(2P), J = 0, 1,2 states are disentangled employing their 
slightly differing invariant masses. The distributions of the event shape observable 
derived from Xft(2P), J = 0, 2 decays or continuum events evolved to the Xfe(2P), J = 1 
scale show clear differences. 

In a similar study [326] events of the type Y(1S) — »• gg7 are compared with 
e^e~ —>■ qq7 events produced in the continuum. The mean charged particle multiplicity 
is measured in both samples as a function of the recoil mass of the photon. The ratio 
of multiplicities in gg and qq systems for recoil masses below 7 GeV is found to be 
compatible with unity. The gluon multiplicities are also shown in figures EEl (a) and (b). 

A method to obtain inclusive gluon jets at high energy has been proposed in [303]. 
The inclusive gluon jet definition is achieved experimentally by selecting hadronic 
decays were both q and q are detected in the same hemisphere of the event. The quarks 
are identified via tagging of heavy quark decays and the gluon jet is identified with the 
other hemisphere of the event. In two studies the mean charged particle multiplicity [328] 
and the charged particle multiplicity distributions [329] were measured in inclusive gluon 
and quark jets. The multiplicity distribution was analysed in terms of higher moments 
and agreement with a NNLO QCD calculation was found for gluon and quark jets. The 
latest result for the ratio of mean charged particle multiplicities in inclusive quark and 
gluon jets of Ej^t = 40.1 GeV is r^h = 1.51 ± 0.02(stat.) ± 0.03(syst.) [330]. This is 
only in qualitative agreement with results from exclusive jets discussed in the previous 
section. 

In [330] many properties of inclusive gluon and quark jets at high jet energies were 
studied. The measurement of the distributions of charged particle energy fractions xe 
allows to test some of the predictions discussed above. From the data for in table 2 
of [330] one can extract {xE)g = 0.044 ± 0.002, {xE)g = 0.0635 ± 0.0003. The ratio is 
found as 

r^^) = {xE)J{xE)g = lAA±0.07 , (64) 

where the error contains an additional contribution of 3% to cover the energy difference 
between the inclusive gluon jets at 40.1 GeV and the quark jets at 45.6 GeV. The 

fratio of 2nd to 0th Fox- Wolfram moment using charged and neutral particles [327] 
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Figure 27. Ratio of scaled momentum distributions of charged particles measured 
in inclusive gluon and quark jets at E-^^t = 40.1 GeV. Superimposed are fits of the 
form (1 — a;)/ ln(l — x) as predicted in [301]. The DELPHI data points arc horizontally 
displaced for clarity. Data from [317,330,331]. 



result for r^^^ is consistent with the result for shown above. The prediction of [301] 
discussed above for the behaviour of the ratio of distributions can be tested by fitting 
a function a(l — x)/ ln(l — x) +6, where a and h free parameters, to the data for xe > 0.2. 
Figure EH presents the results of such fits using the data of [317,330,331]. The fits based 
on the total errors agree with the data; the results are shown in table El The fit to 
OPAL data was performed with b = 0.0 fixed since the fitted value of b was consistent 
with zero. 



Table 14. Results of fits of a(l — a;)/ln(l ~ x) + b as explained in the text to data 
for ratios of FFs from gluon and quark jets. 





a 


b 


xVd.o.f. 


ALEPH yeut = 0.1 


0.176 ±0.064 


0.36 ±0.11 


0.03/1 


DELPHI Merc. 


0.267 ±0.033 


0.067 ±0.057 


4.3/4 


DELPHI Y 


0.306 ±0.010 


0.124 ±0.014 


25/5 


OPAL g.nci 


0.218 ±0.035 


0.0 (fixed) 


1.3/6 



The results for a can be averaged yielding a = 0.295±0.009 with a x^/d.o.f. = 10/3. 
The data quahtatively support the prediction of [301]. However, there are significant 
differences between the fit results indicating substantial systematic differences between 
the different data sets, which have been obtained using different experimental techniques. 
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9. QCD gauge structure 

QCD as the gauge theory of strong interactions assumes that quarks carry one out of 
three strong charges, referred to as colour. The requirement of local gauge symmetry 
under SU(3) transformations in the colour space generates the gauge bosons of QCD: 
an octet of gluons each carrying colour charge and anti-charge, see e.g. [5]. The gluons 
can thus interact with themselves; it is due to this property of the QCD gauge bosons 
that the theory explains confinement and asymptotic freedom via the running of the 
strong coupling as. 

In perturbative QCD at NLO three fundamental vertices contribute: i) the quark- 
gluon vertex with colour factor Cp, ii) the gluon-gluon vertex with colour factor Ca and 
iii) the qq production from a gluon with colour factor TpUf, see e.g. [332]. The colour 
factors are Cp = 4/3, Ca = 3 and TpUj = 1/2 ■ 5 in QCD with SU(3) gauge symmetry 
and specify the relative contribution of the corresponding vertex to observables. In 
NLO QCD the prediction for an observable R is R = Aas + {BcpCp + Bcj^Ca + 
BTpTpnf)Cpa\] NLLA predictions e.g. for event shapes or jet rates decompose in a 
similar way. For NLO predictions for 4-jet observables an analogous decomposition in 
terms of the six possible products of two out of the three colour factors holds [7]; a 
seventh term is negligible and left out. 

Experimental investigations of the gauge structure of QCD are possible because of 
the different angular momenta in the initial and final states of the fundamental vertices. 
It is an important test of QCD to probe the gauge structure in experiments. Several 
techniques with rather different experimental and theoretical uncertainties have been 
developed; we will discuss here some recent results. 

9.1. Four- jet events 

The LEP experiments ALEPH [181] and OPAL [185] have analysed 4-jet final states from 
hadronic Z° decays using the recent QCD NLO predictions (see e.g. [7] and references 
therein). The 4-jet final states are selected by clustering events using the Durham 
algorithm [169] with ?/cut = 0.008 and demanding four jets. At this value of ?/cut the 
4-jet fraction is relatively large (R4 ~ 7%) and the four jets are well separated. 

The energy-ordered 4-momenta pj, i = 1, . . . , 4 of the jets are used to calculate the 
angular correlation observables (see [181,185] for details). As an example, the Bengtsson- 
Zerwas angle xbz is defined by [183] xbz = ^{[Vi x ^2], [Pa x ^4]), i-e. the angle between 
the two planes spanned by the momentum vector pairs (^1,^2) and (^3,^4). Assuming 
that energy ordering selected the primary quarks from the TP decay as (^1,^2) the 
observable xbz is sensitive to the decay of an intermediate gluon to gluons (vertex ii)) 
or quarks (vertex iii)) and the competing process of radiation of a second gluon from a 
primary quark (vertex i)). 

Figure |2H1 (left) presents the uncorrected distribution of cos(xbz) measured by 
ALEPH [181]. Superimposed on the data points is the result of a simultaneous fit of the 
NLO QCD predictions to four angular correlations including xbz and the distribution 
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of the 4-jet rate R4. From the fit values for Q;g(mzo) and the colour factors Ca and Cp 
are extracted. The results from ALEPH are as(mzo) = 0.119 ± 0.027, Ca = 2.93 ± 0.60 
and Cf = 1.35 ± 0.27; from OPAL we have as(mzo) = 0.120 ± 0.023, Ca = 3.02 ± 0.55 
and Cp = 1.34 ± 0.26. The systematic errors are dominated by uncertainties from 
the hadronisation corrections and theoretical errors from estimating missing higher 
order contributions. The hadronisation corrections are implemented using Monte Carlo 
models with standard values for the colour factors. Both analyses use an unconventional 
method of estimating systematic uncertainties. A more conservative evaluation based 
on the information available in [181,185] leads to total errors for ag(mzo), Ca and Cp 
of ±0.048, ±1.06 and ±0.46 for ALEPH and ±0.049, ±1.07 and ±0.47 for OPAL; i.e. 
the errors turn out to be approximately twice as large. 
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Figure 28. (left) Distribution of |cosxbz| before corrections (solid points). 
Superimposed is a fit of NLO QCD (solid lines) [181]. (right) Measurements of 
d£'^(a;E, KH)/dlog(/vH) for gluon (solid points) and quark jets (open circles and solid 
triangles). The lines represent a fit of the LO QCD prediction [333]. 



9.2. Scaling violation in gluon and quark jets 

Jets originating from quarks or gluons should have different properties due to the 
different colour charge carried by quarks or gluons, see section |H| In an analysis by 
DELPHI [333] the scaling violation of the fragmentation function (FF) in gluon and 
quark jets at different energies is compared. From a sample of 3.7 ■ 10^ hadronic 
decays planar 3-jet events with well reconstructed jets are selected using the Durham 
or Cambridge algorithms, but without imposing a fixed value of tjcnt in order to reduce 
biases [315] (see section HH} . In events with the two angles between the most energetic 
jet and the other jets in the range between 100° and 170° a b-tagging procedure is 
applied to the jets. Gluon jets are identified indirectly as the jets without a successful 
b-tag. Jets originating from light (udsc) quarks are taken from events which failed 
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the b-tagging. A correction procedure based on efficiencies and purities determined by 
Monte Carlo simulation yields results for pure gluon or quark jets. 

The jet energy scale is calculated according to kh = -Ejet sin(^/2), where 9 is 
the angle w.r.t. to the closest jet. This definition takes colour coherence effects into 
account [302]. The FF in a jet D^{xe, i^k) is given by the distribution of xe = -Ehadmn/ ^^h 
for the hadrons assigned to the jet. The evolution of the FFs of gluon and quark jets with 
jet energy scale kh is studied in intervals of xe- Figure EH (right) presents the quantity 
dD^{xE, KH)/d log(KH), i.e. the slopes of the scaling violation for a given interval of xe, 
for quark and gluon jets. The steeper slopes corresponding to stronger scaling violations 
of gluon compared to quark jets are clearly visible. A fit of the scaling violations based on 
the LO QCD prediction (DGLAP equation) allows to extract the ratio Ca/Cf resulting 
in Ca/Cf = 2.26 ±0.16. 

9. 3. Event shape fits 

In this analysis [334] the decomposition of the C(ag)+NLLA QCD predictions for event 
shape observables into terms proportional to the colour factors is used. Since the 
sensitivity of event shape distributions measured at LEP 1 alone is not sufficient [335] 
data from a/s = 14 to 189 GeV are used. In this way the colour structure of the 
running of the strong coupling contributes as well. Hadronisation corrections are 
implemented using power corrections, see section 14.5.11 The advantage of using power 
corrections instead of Monte Carlo model based hadronisation corrections is that the 
colour structure of the power corrections is known and can be varied in the fit. 

In the analysis simultaneous fits of as{mzo), Ca and to data for the event shape 
observables 1 — T at ^/s = 14 to 189 GeV and C at ^/s = 35 to 189 GeV are performed. 
The data for 1 — T and C are analysed separately and the results are combined. The 
results are as(mzo) = 0.119 ± 0.010, Ca = 2.84 ± 0.24 and Cf = 1.29 ± 0.18 and 
are shown on figure 1221 below. The errors are dominated by uncertainties from the 
hadronisation correction and from experimental effects. 

9.4- Colour factor averages 

The measurements of the colour factors Ca and Cf or oi x = Ca/Cf discussed above 
or in section IHl can be combined into average values of Ca and Cf taking into account 
correlations between Ca and C^? as well as between different experiments. The variables 
X and y = Tf/Cf are used to define the function 

= Y^i^i- ^)vi'''"\yi -y)+ (65) 

i 

^(Xi - x)v-f{xj -x) + 

Y^ivi - y)^ij\yj - y) , 
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where x and y are the averages, the indices i and j count experiments, vl , vij and 
ffj^ are elements of the inverses of the corresponding covariance matrices for the Xi and 
Ui within experiment i and for the Xj and Ui between experiments. The averages x and 
y are converted to the average values Ca and Cp after the fit is performed. The input 
data for x and y are directly taken from [181,185,320,333] while the results from [334] 
have to be convertedf. 

The covariance matrices between experiments are constructed as follows: the 
ALEPH and OPAL 4-jet analyses have experimental and hadronisation errors partially 
and theory errors fully correlated, the 4-jet analyses and the event shape analysis have 
hadronisation errors partially and theory errors fully correlated, and the 4-jet and event 
shape analyses and the DELPHI FF and the OPAL A^^gg' analyses have their theory 
errors partially correlated. 

The averaging fit is done using only the first term of equation ()65|1 to determine 
the averages Ca and Cp avoiding possible biases from our assumptions on correlations 
of systematic uncertainties between different sets of measurements. The fit is repeated 
using only statistical correlations in the first term to determine the statistical errors and 
using the full covariance matrix to determine the total errors. The final results are 

Ca = 2.89 ± 0.03(stat.) ± 0.21(syst.) (66) 
Cp = 1.30 ± 0.01(stat.) ± 0.09(syst.) 

with correlation coefficient p = 0.82. The relative total uncertainties are about 8% 
for both Ca and Cp. The fit using the full covariance matrix yields Ca = 2.83 and 
Cp = 1.29 in good agreement with our main results. 

Figure 12^1 presents the results of the individual analysis in a. Cp vs. Ca plane 
together with the combined result and the expectations of QCD based on the SU(3) 
gauge symmetry and various other gauge symmetries. The correlation coefficients 
for [181,185] were calculated from the references with the results p = 0.97 (ALEPH) and 
p = 0.93 (OPAL). The error ellipses refer to 90% CL. The combined result is in good 
agreement with the individual analyses and with standard SU(3) QCD while the total 
uncertainties are substantially reduced. The other possibilities for gauge symmetries 
shown on the figure are clearly ruled out. 

10. Conclusions and outlook 

We have shown in this report that hadron production in e"'"e~ annihilation is a useful 
and fruitful environment to test the theory of strong interactions, QCD. The absence of 
interference between initial and final states and the large range of cms energies probed 
by the experiments make many stringent tests of the theory possible. The early results 
from hadron production in e"'"e~ annihilation at low energies were crucial to establish 

fThe results are x = 2.20 ± 0.21(stat.) ± 0.25(syst.), y = 0.388 ± 0.021(stat.) ± 0.051(syst.), 
Pstat. = 0.98 and psyst. = 0.86. 
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Figure 29. Measurements of the colour factors Ca and Cp discussed in this report. 
The ehipses show the correlated measurements using 4-jet events [181,185] or event 
shape distributions [334] while the lines represent the results of determinations of 
Ca/Cf from DELPHI [333] (dashed) and OPAL [320] (solid). The upper sohd and 
dashed lines overlap. The grey filled ellipsis displays the combined result for Ca 
and Cf (see text). The solid triangle and squares show the expectations for various 
assumptions for the gauge symmetry of QCD as indicated on the figure. 

QCD as the theory of strong interactions and thus as an integral part of the standard 
model of high energy physics. 

Studies of differences between quark and gluon jets reveal many properties of the 
gauge bosons of QCD, the gluons, which are correctly predicted by the theory. The 
measurements of jet production rates and event shape observables using theoretically 
and experimentally well behaved observables allow direct tests of advanced perturbative 
QCD predictions and precise determinations of the value of the strong coupling constant 
as- A fundamental prediction by QCD is asymptotic freedom of the coupling at high 
energies and this has been verified directly using data over a large range of cms energies. 
More indirect tests of asymptotic freedom stem from successful comparison of precision 
determinations of as at different energy scales. 

The most reliable and precise determinations of as in e"'"e~ annihilation to hadrons 
employ inclusive observables such as the hadronic branching ratios of the gauge bosons of 
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the electroweak theory, 7, Z° and W"*". The inclusive observables are mostly unaffected 
by hadronisation corrections and they arc accessible to NNLO perturbative QCD 
calculations with 2-loop radiative corrections. The final result for Q;s(mzo) from inclusive 
observables calculated in NNLO perturbative QCD is 

as{mzo) = 0.1211 ± 0.0010(exp.) ± 0.0018(theo.) . 

This result is in good agreement with recent world averages of measurements of the 
strong coupling constant Q;s(mzo) = 0.1184 ± 0.0027 [207] or as{mzo) = 0.1187 ± 
0.0020 [18] which are based on all available analyses. The best measurements of as('^zo) 
including ours shown above yield errors of about 2%. The other determinations of as 
discussed in this report are found to be in good agreement with our main result and we 
note that the results from analyses of the 4-jet rate using the Durham or Cambridge 
algorithms are of similar precision. Further improvements in the near future may be 
expected from NNLO calculations for 3-jet observables which are expected to reduce 
the otherwise dominating theoretical uncertainties and might allow total uncertainties 
of 1% [336]. 

Studies of the gauge structure of QCD using angular correlations in 4-jet final 
states, 3-jet observables at different energy scales and quark and gluon jet differences 
were discTisscd. The results for the colour factors Ca and Cp were observed to be in 
good agreement and combined values were calculated. The results are 

Ca = 2.89 ± 0.03(stat.) ± 0.21(syst.) 
Cf = 1.30 ± 0.01(stat.) ± 0.09(syst.) 

with correlation p — 0.82. These results are in good agreement with expectations from 
QCD with the SU(3) gauge symmetry = 3 and Cp = 4:/3. 

In summary, e"'"e~ annihilation to hadrons has been and still is indispensable to 
stringently probe the theory of strong interactions. QCD has been established as a 
part of the standard model of high energy physics and has been tested to a precision of 
2% with some well suited observables. This impressive result has been achieved by the 
LEP experiments, SLD, the re-analysis of JADE data and the parallel theoretical work 
stimulated by the availability of precise data. Future high energy physics programmes 
such as the LHC or the ILC will benefit greatly from this achievement, because precise 
and reliable QCD predictions will be essential for disentangling possible new physics 
signals from standard processes. 
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